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CROSS-REFERENCE TO RELATED APPLICATION 

This application claims priority under 35 U.S.C. § 1 19(e) to Provisional 
U.S. Patent Application Serial No. 60/413,205, filed September 23, 2002, the entire 
contents of which are incorporated herein by reference. 

FIELD OF THE INVENTION 

The present invention relates to estimating precise position of a stationary or 
moving object using multiple satellite signals and a network of multiple receivers. The 
present invention is particularly suited to position estimation in real-time kinetic 
environments where it is desirable to take into account the spatial distribution of the 
ionosphere delay. 

BACKGROUND OF THE INVENTION 

Satellite navigation systems, such as GPS (USA) and GLONASS (Russia), are 
intended for accuracy self-positioning of different users possessing special navigation 
receivers. A navigation receiver receives and processes radio signals broadcast by 
satellites located within line-of-sight distance, and from this, computes the position of 
the receiver within a pre-defined coordinate system. However, for military reasons, the 
most accurate parts of these satellite signals are encrypted with codes only known to 
military users. Civilian users cannot access the most accurate parts of the satellite 
signals, which makes it difficult for civilian users to achieve accurate results. In 
addition, there are sources of noise and error that degrade the accuracy of the satellite 
signals, and consequently reduce the accuracy of computed values of position. Such 
sources include carrier ambiguities, receiver time offsets, and atmospheric effects on the 
satellite signals. 

The present invention is directed to increasing the accuracy of estimating the 
position of a rover station in view of carrier ambiguities, receiver time offsets, and 
atmospheric effects. 



SUMMARY OF THE INVENTION 

The present invention relates to a method and apparatus for determining the 
position of a receiver (e.g., rover) with respect to the positions of at least two other 
receivers (e.g., base receivers) which are located at known positions. The knowledge of 
the precise locations of the at least two other receivers (located at known positions within 
the coordinate system) makes it possible to better account for one or all of carrier 
ambiguities, receiver time offsets, and atmospheric effects encountered by the rover 
receiver, and to thereby increase the accuracy of the estimated receiver position of the 
rover (e.g., rover position). 

In a first aspect of the present invention, an exemplary apparatus/method 
comprises receiving the known locations of a first base station and a second base station, 
obtaining a time offset representative of the time difference between the clocks of the 
first and second base stations, and obtaining measured satellite data as received by the 
rover, the first base station, and the second base station. The measured satellite data 
comprises pseudo-range information. The exemplary apparatus/method generates a first 
set of residuals of differential navigation equations associated with a set of measured 
pseudo-ranges related to a first baseline (R-Bl) between the rover and the first base 
station. The residuals are related to the measured satellite data received by the rover 
station and the first base station, the locations of the satellites, and the locations of the 
rover station and the first base station. The exemplary apparatus/method also generates a 
second set of residuals of differential navigation equations associated with a set of 
measure pseudo-ranges related to a second baseline (R-B2) between the rover and the 
second base station. These residuals are related to the measured satellite data received 
by the rover station and the second base station, the locations of the satellites, and the 
locations of the rover station and the second base station. The exemplary 
apparatus/method estimates the rover's location from the first set of residuals, the second 
set of residuals, the time offset between the clocks of the first and second base stations. 

In a second aspect of the present invention, the position information and 
measured satellite data from a third base station are received. Also, a time offset 
representative of the time difference between the clocks of the first and third base 
stations is obtained. Thereafter, the exemplary apparatus/method also generates a third 
set of residuals of differential navigation equations associated with a set of measured 
pseudo-ranges related to a third baseline (R-B3) between the rover and the third base 
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station, and estimates the rover's location further with the third set of residuals and the 
time offset between the clocks of the first and third base stations. 

In a third aspect of the present invention, which may be applied with any of the 
aspects of the present invention described herein, the exemplary apparatus/method 
generates the time offset(s) for the base stations from the measured satellite data of the 
base stations provided to it. 

In a fourth aspect of the present invention, which may be applied with any of the 
aspects of the present invention described herein, the exemplary apparatus/method 
obtains measured satellite carrier phase data as received by the rover and the first base 
station (the first base line). The exemplary apparatus/method generates a fourth set of 
residuals of differential navigation equations associated with a set of measured carrier 
phase data related to the first base line, and resolves the cycle ambiguities in the carrier 
phase data from the fourth residual and one or more of the first, second, and third 
residuals. The resolved cycle ambiguities may take the form of floating ambiguities, 
fixed-integer ambiguities, and/or integer ambiguities. The exemplary apparatus/method 
estimates the rover's location further with the fourth set of residuals and the resolved 
cycle ambiguities associated with the first base line. 

In a fifth aspect of the present invention, which may be applied with any of the 
aspects of the present invention described herein, the exemplary apparatus/method 
obtains measured satellite carrier phase data as received by the second base station, and 
further obtains a set of cycle ambiguities related to a set of satellite phase measurements 
associated with the baseline between the first and second base stations (B1-B2). The 
exemplary apparatus/method generates a fifth set of residuals of differential navigation 
equations associated with a set of measured carrier phase data related to the second base 
line between the rover and the second base station (R-B2), and resolves the cycle 
ambiguities in the carrier phase data from the fifth residual, the set of cycle ambiguities 
associated with the baseline between the first and second base stations, and one or more 
of the first, second, third, and fourth residuals. The resolved cycle ambiguities may take 
the form of floating ambiguities, fixed-integer ambiguities, and/or integer ambiguities. 
The exemplary apparatus/method estimates the rover's location further with the fifth set 
of residuals and the resolved cycle ambiguities associated with the second base line (R- 
B2). 

In a sixth aspect of the present invention, which may be applied with any of the 
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aspects of the present invention described herein, the exemplary apparatus/method 
receives measured satellite carrier phase data as received by the third base station, and 
further obtains a set of cycle ambiguities related to a set of satellite phase measurements 
associated with the baseline between the first and third base stations (B1-B3). The 
exemplary apparatus/method generates a sixth set of residuals of differential navigation 
equations associated with a set of measured carrier phase data related to the third base 
line between the rover and the third base station (R-B3), and resolves the cycle 
ambiguities in the carrier phase data from the sixth residual, the set of cycle ambiguities 
associated with the baseline between the first and third base stations, and one or more of 
the first, second, third, fourth, and fifth residuals. The resolved cycle ambiguities may 
take the form of floating ambiguities, fixed-integer ambiguities, and/or integer 
ambiguities. The exemplary apparatus/method estimates the rover's location further with 
the sixth set of residuals and the resolved cycle ambiguities associated with the third base 
line (R-B3). 

In a seventh aspect of the present invention, which may be applied with any of 
the aspects of the present invention described herein, the exemplary apparatus/method 
obtains a first set of first ionosphere delay differentials associated with the satellite 
signals received along the base line formed by the first and second base stations, and 
generates ionosphere delay corrections to one or more of the above-described first 
through sixth residuals from the first set of first ionosphere delay differentials, the 
locations of the first and second base stations, and an estimated location of the rover 
station. As a further option, the exemplary apparatus/method obtains set of second 
ionosphere delay differentials associated with the satellite signals received along the base 
line formed by the first and third base stations (or the base line formed by the second and 
third base stations), and generates the ionosphere delay corrections to the one or more of 
the above-described first through sixth residuals further from the second set of 
ionosphere delay differentials and the location of the third base station. 

In an eight aspect of the present invention, which may be applied to the seventh 

and further aspects of the present invention described herein, the exemplary 

apparatus/method forms one or more of the first through sixth residuals to account for 

second order effects in the ionosphere delay corrections applied to the baselines 

associated with the rover, and further generates an estimate of the second order effects. 

In a ninth aspect of the present invention, which may be applied with any of the 
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aspects of the present invention described herein, the exemplary apparatus/method 
generates the time offset representative of the time difference between the clocks of the 
first and second base stations, and optionally the time offset representative of the time 
difference between the clocks of the first and third base stations. As a further option, the 
5 apparatus/method generates the time offset representative of the time difference between 
the clocks of the second and third base stations, and performs a consistency check of the 
time offsets. 

In an tenth aspect of the present invention, which may be applied with several of 
the aspects of the present invention described herein, the exemplary apparatus/method 

10 generates the resolved cycle ambiguities associated with the baseline between the first 
and second base stations (B1-B2), and optionally generates the resolved cycle 
ambiguities associated with the baseline between the first and third base stations (BI- 
BS). As a further option, the apparatus/method generates the resolved cycle ambiguities 
associated with the baseline between the second and third base stations (B2-B3), and 

15 performs a consistency check of the cycle ambiguities associated with the baselines 
between the base stations. 

In an eleventh aspect of the present invention, which may be applied with any 
aspects of the present invention described herein which account for ionosphere delays, 
the exemplary apparatus/method generates the first set of first ionosphere delay 

20 differentials, and optionally the second set of first ionosphere delay differentials. As a 
further option, the apparatus/method generates a third set of ionosphere delay 
differentials associated with the satellite signals received along the base line formed by 
the second and third base stations, and generates therefrom three sets of ionosphere delay 
differentials which are self consistent. 

25 Accordingly, it is an object of the present invention to increase the accuracy of 

estimating the rover's position using information from one or more baselines associated 
with two or more base stations. 

It is a further object of the present invention to enable the spacing between the 
rover and the base stations to increase while maintaining or improving the accuracy of 

30 the estimation of the rover's position. 

These and other objects of the present invention will become apparent to those 
skilled in the art from the following detailed description of the invention, the 
accompanying drawings, and the appended claims. 



BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a perspective view of a rover station (R) and three base stations (Bl, B2, 
B3) in an exemplary network according to exemplary embodiments of the present 
invention. 

FIG. 2 is a top-plan schematic drawing the rover station (R) and three base 
stations (Bl, B2, B3) in the exemplary network shown in FIG. 1 according to the present 
invention. 

FIG. 3 is a perspective view of the ionosphere delay differentials to selected 
exemplary methods according to the present invention. 

FIG. 4 is a top-plan schematic view of a road application where a reduced level of 
interpolation of ionosphere delays may be used according to the present invention. 

FIG. 5 is a schematic diagram of an exemplary rover station according to the 
present invention. 

FIG. 6 is a general flow diagram of embodiments of the present invention. 
FIG. 7 is a schematic diagram of an exemplary computer program product 
according to the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

FIG. 1 is a perspective view of a rover station (R) and three base stations (Bl, B2, 

B3) in an exemplary network according exemplary embodiments of the present 

invention, and FIG. 2 is a top-plan schematic drawing thereof. The present invention 

pertains to estimating the position of the rover station with information provided from 

two or three of the base stations and with satellite measurements made by the rover 

station. Referring to FIG. 1, each station has a receiver that receives the satellite 

positioning signals with a satellite antenna (shown as a substantially flat disk). A 

plurality of satellites S1-S4 are depicted in FIG.l, with the range between each satellite 

and each antenna being depicted by a respective dashed line. In the example shown and 

FIG. 1, the Rover station is operated by a human user, and has a positioning pole for 

positioning the Rover's satellite antenna over a location whose coordinates are to be 

determined. The Rover's satellite antenna is coupled to the receiver's processor, which 

may be disposed in the backpack and carried by the user. The user interacts with the 

receiver's processor through a keypad/display. The receiver also has a radio modem 

(more formally a demodulator) that can receive data from the base stations through a 
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conventional RF antenna and relay it to the processor. In other embodiments, the 
Rover's satellite antenna may be mounted to a vehicle, and the receiver may operate 
independently and automatically without the need of a human operator. The present 
invention is applicable to these and other physical embodiments. 
5 The satellite signals comprise carrier signals which are modulated by pseudo- 

random binary codes, which are then used to measure the delay relative to local reference 
clock or oscillator. These measurements enable one to determine the so-called pseudo- 
ranges between the receiver and the satellites. The pseudo-ranges are different from true 
ranges (distances) between the receiver and the satellites due to variations in the time 

10 scales of the satellites and receiver and various noise sources. To produce these time 

scales, each satellite has its own on-board atomic clock, and the receiver has its own on- 
board clock, which usually comprises a quartz crystal. If the number of satellites is large 
enough (four or more), then the measured pseudo-ranges can be processed to determine 
the user location (e.g., X, Y, and Z coordinates) and to reconcile the variations in the 

1 5 time scales. Finding the user location by this process is often referred to as solving a 
navigational problem or task. 

More specifically, the GPS system employs a constellation of satellites in orbit 
around the earth at an altitude of approximately 26,000 km. Each GPS satellite transmits 
microwave radio signals in two frequency bands located around 1575.42 MHz and 

20 1227.6 MHz, referred to as LI band and L2 band, respectively. The GPS LI -band signal is 
modulated by a coarse/acquisition code (C/A) and a precision ranging code (P-code). The 
L2-band signal is binary modulated by the P-code. The GPS C/A code is a pseudo-random 
(PR) Gold code that is specific to each satellite and is used to identify the source of a 
received signal. The P-code is a pseudo-random code signals and is also specific to each 

25 satellite, having a symbol rate which is ten time more than C/A, which reduces the 

granularity by a factor of ten. The GPS satellite transmission standards are set in detail by 
the ICD-GPS-200, Revision C, ARINC Research Corporation, 10 October, 1993. 

The satellites of the GLONASS system transmit signals in two frequency bands 
around 1602 MHz, and around 1246 MHZ, referred to also as LI band and L2 band, 

30 respectively. The GLONASS LI -band signal is modulated by a pseudo-random C/A code 
and a pseudo-random P-code. The L2-band signal is modulated by the P-code. Unlike 
GPS, in which all of the satellites transmit on the same nominal frequency, the GLONASS 
satellites each transmit at a unique frequency in order to differentiate between the satellites. 



The GLONASS Ll-band carrier frequency is equal to 1602 MHz + /*0.5625 MHz, where / 
is an integer index ranging between 0 and 24 that identifies the satellites. The GLONASS 
L2-band carrier frequency is equal to 1246 MHz + 7*0.4375 MHz. Details of the 
GLONASS signals maybe found in GLONASS ICD, Version 4.0 f 1998. 

The distance between a receiver and a satellite (called the "receiver-to-satellite 
range*') is determined by measuring the time that it takes for the signal to pass from the 
satellite to the receiver, provided that the position of the satellite is known. The satellites 
and receivers have internal clocks that are synchronized to a single GPS time. For each 
satellite signal being tracked, the receiver generates a local version of the satellite's 
expected PR-code, and then retards that version in time until the local version correlates 
{i.e., matches) with the received satellite signal. Thereafter, the satellite signal is tracked by 
advancing or retarding the local version of the PR-code. The carrier phase of the satellite 
can also be tracked, which is usually done by tracking the Doppler shift of the satellite 
signal. The positions of the satellites are, except for minor variations, highly predictable as 
a function of time, and the receiver generally carriers a model of the satellite's position as a 
function of GPS time. In theory, by determining the ranges to three different satellites, the 
receiver can perform a three-dimensional triangulation to find its position. But because 
limitations in the accuracy of the receiver's clock, the internally generated time is offset 
somewhat from true GPS time. Thus, the ranges to at least four different satellites are 
simultaneously measured in order to be able to solve for four unknowns, namely the three 
coordinates of the position of the receiver location {e.g., x, y, and z) and an offset of the 
receiver clock time from the GPS time. The location is usually performed with respect the 
defined Cartesian coordinates frame. 

In theory, a GPS receiver can track both the C/A code and the P-code of a satellite. 
The C/A code generally provides accuracy to within 20 - 50 meters, and the P-code 
generally provides greater accuracy to within 10 meters because of its higher repetition rate 
(less granularity) relative to the C/A code. However, knowledge of the P-code is restricted 
to military users and not available to civilian users. Nonetheless, some approaches for 
tracking P-codes have been developed. 

The desire to guarantee the solution of navigational tasks with accuracy better than 

10 meters, and the desire to raise the stability and reliability of measurements, have led 

to the development of the mode of "differential navigation ranging," also called 

"differential navigation" (DN). In the DN mode, the task of finding the user position is 
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performed relative to a Base station (Base), the coordinates of which are known with the 
high accuracy and precision. The Base station has a navigation receiver that receives the 
signals of the satellites and processes them to generate measurements. The results of 
these measurements enable one to calculate corrections, which are then transmitted to a 
roving GPS receiver, which the user has set up. We call this GPS receiver the "Rover 
station," or "Rover receiver." By using these corrections, the roving GPS receiver gains 
the ability to compensate for the major part of the strongly correlated errors in the 
measured pseudo-ranges, and to substantially improve the accuracy of the estimate of its 
position. 

Usually, the Base station is immobile during measurements. The rover station may 
be either immobile or mobile. Depending on the navigational tasks to be solved, 
different modes of operation may be used in the DN mode. They differ in the way in 
which the measurement results are transmitted from the Base to the Rover. In the post- 
processing (PP) mode, these results are transmitted as digital recordings and go to the 
user after all the measurements have been finished. In the PP mode, the user reconstructs 
his or her location for definite moments in the past. 

Another mode is the Real-Time Processing (RTP) mode, and it provides for the 
positioning of the Rover receiver just during the measurements. The RTP mode uses a 
communication link (such as the radio communication links shown in FIG. 1), through 
which all the necessary information is transmitted from the Base to the Rover receiver in 
digital form. 

Further improvement of accuracy of differential navigation may be reached by 
supplementing the measurements of the pseudoranges with the measurements of the 
phases of the satellite carrier signals. If one measures the carrier phase of the signal 
received from a satellite in the Base receiver and compares it with the carrier phase of the 
same satellite measured in the Rover receiver, one can obtain measurement accuracy to 
within several percent of the carrier's wavelength, i.e., to within several centimeters. 

The practical implementation of those advantages, which might be guaranteed by 
the measurement of the carrier phases, runs into the problem of there being ambiguities 
in the phase measurements. 

The ambiguities are caused by two factors. First, the difference of distances AD 
from any satellite to the Base and Rover is usually much greater than the carrier's 

wavelength A. Therefore, the difference in the phase delays of a carrier signal A<p=AD/A 
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received by the Base and Rover receivers may substantially exceed one cycle. Second, it 
is not possible to measure the integer number of cycles in A<p from the incoming satellite 
signals; one can only measure the fractional part of A<p. Therefore, it is necessary to 
determine the integer part of A<p, which is called the "ambiguity". More precisely, we 
need to determine the set of all such integer parts for all the satellites being tracked, one 
integer part for each satellite. One has to determine this set along with other unknown 
values, which include the Rover's coordinates and the variations in the time scales. 

In a general way, the task of generating highly-accurate navigation measurements 
may be formulated by defining a set of unknowns and system of relationships between 
the unknowns and measured navigation parameters. The vector of unknowns, denoted 
herein as include three Rover coordinates (usually along Cartesian axes X, Y, Z) in a 
given coordinate system (sometimes time derivatives of coordinates are added too); the 
variations of the time scales which is caused by the phase drift of the local main 
reference oscillator; and n integer unknown values associated with the ambiguities of the 
phase measurements of the carrier frequencies. The value of n is determined by the 
number of different carrier signals being processed, and accordingly coincides with the 
number of satellite channels actively functioning in the receiver. At least one satellite 
channel is used for each satellite whose broadcast signals are being received and 
processed by the receiver. Some satellites broadcast more than one code-modulated 
carrier signal, such as a GPS satellite, which broadcasts a carrier in the Li frequency 
band and a carrier in the L 2 frequency band. If the receiver processes the carrier signals 
in both of the Li and L 2 bands, the number of satellite channels (n) increases 
correspondingly. 

Two sets of navigation parameters are measured by the Base and Rover receivers, 
respectively, and are used to determine the set of unknowns, which is sometimes called 
the "state vector." Each set of navigation parameters includes the pseudo-range of each 
satellite to the receiver, and the fUll (complete) phase of each satellite carrier signal. 
Each pseudo-range is obtained by measuring the time delay of a code modulation signal 
of the corresponding satellite (C/A code or P-code). The code modulation signal is 
tracked by a delay-lock loop (DLL) circuit in each satellite-tracking channel. The full 
phase of a satellite's carrier signal is tracked by a phase-lock-loop (PLL) in the 
corresponding satellite tracking channel. (The DLL and PLL are, for example, provided 
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by the demodulator 120 of an exemplary rover shown in FIG. 5). An observation vector 
is generated as the collection of the measured navigation parameters for specific 
(definite) moments of time. 

The relationship between the state vector and the observation vector is defined by a 
well-known system of navigation equations. Given an observation vector, the system of 
equations may be solved to find the state vector if the number of equations equals or 
exceeds the number of unknowns in the state vector. In the latter case, conventional 
statistical methods are used to solve the system: the least squares method, the method of 
dynamic Kalman filtering, the method of least squares and various modifications of these 
methods. 

One method may comprise the following steps. The measured values of the 
pseudo-ranges and full phases at specific (definite) moments of time, along with an 
indication of the satellites to which these measurements belong to and the time moments 
of the measurements, are transmitted from the Base to the Rover (such as through the 
communication link or as recordings). Corresponding values are measured in the Rover 
receiver. The processing includes the determination of the single differences of the 
pseudo-ranges and full phases between the Base and Rover measurements for each 
satellite. The strongly correlated errors are largely compensated (i.e., substantially 
cancelled) in the single differences. Then, the residuals of the single differences are 
calculated by subtraction of calculated values from the measured results. The processing 
of residuals allows one to linearize the initial system of navigation equations (sometimes 
several subsequent iterations are necessary for that), which makes possible the use of the 
well-developed body of mathematics for solving systems of linear equations. The 
subsequent iterative solution of the linearized system of navigation equations is 
equivalent to the minimization of the sum of squared non-linear residuals by the Gauss- 
Newton minimization method. The components of the state vector, with the n 
ambiguities included, are found as a result of the solution. But the calculated values of 
the ambiguities are not necessarily integer numbers, and are often floating point 
numbers. Because of this, they are called "float ambiguities," or "floating ambiguities," 
at this stage of the solution. To find true values of the integer ambiguities one uses the 
procedure of rounding off the float ambiguity vector to the nearest set of integers. This 
process is called 'the ambiguity resolution'. After the ambiguity resolution has been 

done, is it possible to generate more accurate values of residuals and then, by solving the 
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system of equation again, to find the coordinate values for the baseline connecting the 
Base and Rover, and consequently to more accurately estimate the coordinates of Rover 
and the correction to its clock drift. 

There are fractional parts associated with each ambiguity due to the time offset in the 
5 receiver's carrier generator, and as a practical matter, it is more convenient to include this 
fractional offset into the floating ambiguities and "integer" ambiguities. In this case, one 
additionally solves for the unknown fractional offset of the receiver's carrier generator, and 
then resolves the "integer" ambiguities in the form of integers plus a common fractional 
part that is related to the time offset of the receiver's carrier generator. We call these 

1 0 resolved ambiguities "fixed-integer ambiguities" because, although these ambiguities have 
fractional parts, the difference between any two fixed-integer ambiguities for two 
corresponding satellites measured by the same pair of receivers is an integer. 

Even with all of the-above described processing to account for clock offsets and 
carrier ambiguities, there are additional factors which affect the accuracy of measurements 

15 made with GPS and/or GLONASS signals. As one factor, the trajectory of each satellite (or 
its initial data called the ephemeris), is elliptical and is affected by natural causes such as 
solar winds. The accuracy of any measurement is dependent upon knowledge of the 
position of the satellites at certain time. An estimate of the ephemeris is calculated on earth 
for each of the satellites and is periodically uploaded to the satellite. The position 

20 information of a satellite is encoded onto a low frequency (50 Hz) signal which is 

modulated on to one of the carrier signals, and transmitted to the GPS receiver on earth. 

Two additional and important factors that affect the accuracy of measurements 
made with GPS and/or GLONASS signals are the effects of the troposphere and ionosphere 
on propagation of signals from the satellites to the receivers. The troposphere is the lower 

25 part of the atmosphere and variations in the temperature, pressure, and humidity lead to 
spatial variations in the signal propagation. The ionosphere is at the upper part of the 
atmosphere and has a slice of ionized gas at the altitude of around 300 km. The density of 
ionized particles is sufficiently high to affect the propagation of electromagnetic signals, 
and has a spatial variation and a time variation. The ionosphere effect becomes even more 

30 important during years of high solar activity. These variations in the troposphere and 

ionosphere introduce variable delays in the propagation the satellite signals to the receivers. 

If the base and rover are widely separated, these delays will be significantly different for the 

base and rover stations, and will introduce error into the estimation of the rover's position. 
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With that general overview, we now describe the invention in greater detail. For 
the sake of simplifying the presentation, and without loss of generality, we assume 'W 
satellites identified by the index "s", s = 1, . . N 9 and we assume that each receiver can 
track the satellite's LI band signal and the satellites L2 band signal. The receivers 
extract timing information from the satellite signals, and report this information a 
predefined increments k of time, which we call epochs k. The time between epochs can 
be selected by the user, and generally ranges between 0.1 seconds to 2 seconds, with 1 
second being typical. The clocks of all of the receivers are typically accurate to within 
several milliseconds of the true GPS time, and for practical purposes the receivers can 
determined the number k of the current epoch from their clocks. The following timing 
information can be extracted from each satellite V at each epoch by each receiver 

1 . The low- frequency (50 Hz) information signal which provides information on 
the orbit (position) of the satellites which enables the receiver to determine 
the satellite's position. 

2. the pseudo-range pfjf derived from the LI -band C/A code (or optionally the 

LI -band P-code for military users), and having the units of length (e.g., 
meters [m]); 

3. the pseudo-range p^ s derived from the L2-band C/A code (or optionally the 

L2-band P-code for military users), and having the units of length (e.g., 
meters [m]); 

4. the LI -band carrier phase <p^ s derived from the LI -band carrier, and 

typically having the dimensionless units of cycles (but is sometimes 
expressed in terms of meters by multiplying the cycles by the LI -band 
wavelength X L1 ' S ; and 

5. the L2-band carrier phase s derived from the L2-band carrier, and 

typically having the dimensionless units of cycles (but is sometimes 
expressed in terms of meters by multiplying the cycles by the L2-band 
wavelength X L2,S . 
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Each of pQ s , pf? 2 k > s , and </>^ s are known as observed quantities, or 

"observables," because the are measured by the receiver from the satellite signals. The 
observables are related to the speed of light "c," the wavelengths X L1,S and >P> S of the Ll- 
band and L2-band carriers of satellite "s the frequencies/ 1 ^ 1 ' 5 =c/ X L1 ' S and/ 1 - 2 ' 5 ^/X 12 '* 
of the LI -band and L2-band carriers of satellite "5," various variables that are to be 
determined, and various noise sources, in the following manner: 

Pr)k S = K,k -*r,k) + T? >k + + nft [1A] 



\ A J 
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where variables (which are to be determined) are: 

• R^k is the true range between the satellite V and receiver V ' at the epoch "it", 
having dimension of distance (usually expressed as meters [m]); 

• r/ is the time offset from true GPS time of the clock of satellite V at epoch 

• T r ,k is the time offset from true GPS time of the clock of receiver V at epoch "A"; 

• T* k is the delay due to tropospheric effects between satellite "s" and receiver "r" 
at epoch 

• I s r k is the delay due to ionospheric effects between satellite "s" and receiver "r" 
at epoch (the full amount of I s r k affects the pseudo-range in the LI -band); 

• Ny X * s is the integer phase ambiguity in the LI -band observable <j> L \ s ; 

• Nf? 2 ' s is the integer phase ambiguity in the L2-band observable <f> L \ s ; 
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• if/, 1 is the initial phase offset of the LI -band phase tracker of receiver "r"; and 

• \/fr 2 is the initial phase offset of the L2-band phase tracker of receiver "r", 



and where various noise sources are: 

• n r,k 1S * e total e ff ect of the noise sources on the pseudo-range measurement at 

epoch in receiver V when tracking the LI -band C/A-code signal (or P-code 
signal) of satellite 'V; 

• n rjc' S * s * e tota ' ef f ect of the noise sources on the pseudo-range measurement at 

epoch in receiver V when tracking the L2-band C/A-code signal (or P-code 
signal) of satellite 'V; 

• v^g 5 is the total effect of the noise sources on the carrier-phase measurement at 
epoch "A" in receiver V when tracking the LI -band carrier of satellite V; and 

• v^jf is the total effect of the noise sources on the carrier-phase measurement at 

epoch "A" in receiver V when tracking the L2-band carrier of satellite V. 
Forms [1A]-[1D] are known as the satellite navigation equations. 

In embodiments of the present invention, we will be forming differences between 
the observables of pairs of stations, and differences between instances of the satellite 
navigation equations. The vector, or straight line, between a pair of two such stations is 
called a baseline. The baselines between all of the stations shown in FIG. 1 are best seen 
in the top plan view of FIG. 2. We generally note the difference of an observable 
measured by two stations with the prefix notation A q r , where each of "q" and "r" 

represents one of the stations Bl, B2, B3, and R. For example, 

A 1,rP k -P r ^ ~P q ; k [2A] 

^Pp'-Pfr-Pft [2B] 

**Ai %M -9%f [2D] 

Forms [2 A] -[2D] are often referred to as the between-station (or between-receiver) single 
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differences of the signals of satellite "s" between stations (receivers) V and "q n m The 
between station operator A gr can be applied to the variables and noise sources as well. 

For example, A q , r R* k = R s r>k and A q>r n£ l >* =n% s -nffi. With this 

background, we can relate the between-station single differences of the observables to 
the differences in the variables as follows: 

* g ,rP£ hs = -cA q>r r k + A g , r T£ + A,, r /* + A q , r n^ s [3 A] 



A 4,rP£ 2 ' S = A q ,rK -cA q , r T k + A ff>r lJ + 



( ^ s 

7 L\,s 



2 



A q>r I s k+ A q>r nj; 2 > s [3B] 



Vtf' 5 = -jnjWi -f LX ' SA q ,m +-^-A f . r i)f 



-ij A,, r /f + A qr N Ll ' s + A qr ijf Li + A^vf * 



[3C] 



[3D] 



10 The difference form of [3 A] is generated by forming two instance of form [1 A] for 

receivers V and * y , and then subtracting the two instances (the instance for receiver 
'V is subtract from the instance for receiver V). Difference forms [3B], [3C], and 
[3D], are formed in a similar manner from corresponding instances of forms [IB], [1C], 
and [ID], respectively. Forms [3A]-[3D] are the single-differences of the navigation 

15 equations. The benefit of forming the between-station single differences is that the error 
term representing the time offset of the satellite clock, r/, is cancelled out in the 
differences. We emphasize that the forms [3 A]-[3D] can be applied to each satellite that 
can be observed by receivers "q n and "r" Higher-order differences of the navigation 
equations, such as the double-differences of the navigation equations, can be formed and 

20 are known to the art. For example, a common double-difference equation is the 

difference between two singles-difference equations associated with a common baseline, 
but with each single-difference equation being based on a different satellite. The present 
invention may also use these higher-order differences of the navigation equations, 
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although the single differences of the navigation equations are currently preferred. Each 
of these difference forms is generically referred to as a differential navigation equation. 

Between Base Station Processing, Part I. 

5 Because the locations of the base stations are known, and because the locations of 

the satellites as a function of epoch time (k) are also known, the between station 

differences of each true range A qr R s k = R s rk -R s qk associated with each satellite 'V 

can be generated from existing information. If vectors X r , X q , and X s represent the 
locations of receiver V\ receiver "q 9 \ and satellite V, respectively, in a coordinate 
10 system (e.g., Cartesian coordinate system), and the operator || ■ || represent the distance 

operator in that coordinate system, then A qr R s k = || X s - X r || - || X s - X q ||. This reduces 

by one the number of unknowns in forms [3A]-[3D]. In addition, the Goad-Goodman 
model may be used to model the difference in troposphere effects (with an error less than 

a few percent), and thus the difference A q r T^ may be estimated based on the positions 

15 of the receivers and the satellite (more specifically, the angles between the receivers and 

satellites). The noise sources A qr n Ll ' s , A qr n L2 > s , A qr v LX * s , and A q r v L2 > s can 

never be known, but they are generally zero-mean and their effect can be reduced by 
averaging. Thus, the number of solvable unknowns in forms [3 A] -[3D] may be reduced 

to the following six: (1) A^, (2) A^ r /|, (3) A^N 11 ** , (4) A^ L1 , (5) 

20 A q r N L2iS (6) A q r \j/ L1 . The first unknown varies with time and is common to all of 

the satellites being tracked by receivers "q" and "r." The second unknown varies with 
time and is specific to the satellite "s" being tracked by the pair of receivers "q n and "r." 
The third and fifth unknowns are each specific to the satellite V being tracked by the 
pair of receivers "q n and "r", and each does not normally vary with time unless a cycle 
25 slip occurs in the receivers phase-lock loop. The forth and sixth unknowns are specific 
to the pair of receivers "q" and "r," and do not normally vary with time. 

As a practical matter, we note that the third and fourth unknowns are related to 
the LI -carrier phase measurement of satellite ' V and are time-independent (unless a 
cycle slip occurs). In theory, these unknowns can be represented as a combined 
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unknown A qr N Ll ' s = A qr N Ll > s + A qr y/ Ll , which we call the between-station "fixed- 
integer ambiguity" for the LI carrier of satellite "s," as measured at stations "q n and "r." 
We call these resolved ambiguities "fixed-integer ambiguities'' because, although these 
ambiguities have fractional parts, the difference between any two fixed-integer 
ambiguities for two corresponding satellites measured by the same pair of receivers is an 
integer. In a similar manner, the fifth and sixth unknowns are related to the L2-carrier 
phase measurement of satellite V and are time-independent (unless a cycle slip occurs); 

they can be represented as a combined unknown A n r N L2 * s = A n r N L2 > s + A„ rW 12 * 

which we call the between-station "fixed-integer ambiguity" for the L2 carrier of satellite 
"s", as measured at stations "q" and "r." 

As a further practical matter, one first generates floating point versions of the 
fixed-integer ambiguities, and then applies a discretization process to the floating point 
version to generate the fixed-integer ambiguities. We call these versions the floating- 
point ambiguities, and denote them as A qr N Ll > s s A qr N Ll > s + A q r y/ LX and 

Thus, the solvable unknowns for each instance of form [3A]-[3D] can be reduced 
to four: (1) A q9r t k9 (2) A^/f, (3) A q ^ r N Ll > s {A q , r N Ll > s ), and (4) A qr N L2 ' s 

( A q r N L2 ' s ). The first unknown is common to all the satellites being tracked by the 

pair of receivers "q" and "r," whereas the last three are specific to the satellite "j." The 
first two unknowns vary with time, whereas the last two are time invariant (unless a 
cycle slip occurs). For the base stations, the known values of forms [3 A] - [3D] can be 



represented as: 








A h LX,s 


= A n Ll > s - 


~ A q,r R k ~ A q,r T k 


[4A] 


A h L2 > s 


- a n L2,s 

^q.rPk 




[4B] 


A n Ll > s 


= A q,rh 


1 A nS 1 a qrS 

A Ll,s ^q> r *k xL \ yS a q,r J k 


[4C] 


A n L2 > s 




x L2,s A 1> rR k x L2,s ^1> rT k 


[4D] 



Forms [4A]-[4D] comprise the known terms of the single-differences of the navigation 
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equations, and we call them the "residuals" of the single differences, or more generally 
the "residuals of a set of differential navigation equations". The magnitude of each 
residual is generally less than the magnitudes of one or more of the known terms that 
forms the residual. 

Then, the unknowns can be estimated from the residuals using the following 
forms by any of number of solution methods known to the art: 



* q ,rPi l ' S =-f LX > S ± q ,rT k ~ A q ,rI S k * \, r N^ + A q ^ 



u 



L2,s 



Ll,s 



,L\,s 



[5A] 
[5B] 
[5C] 



a n L2 ' s --f L2 > s A t, 



,L2,s 



A q , r q + A qtr N L2 ' s + A q , r v^>* [5D] 



L2,s 



All of the residuals, forms, known values, and unknown variables are real-valued. 
The estimation process usually includes instances of forms [5A]-[5D] for several 
satellites, and for several epochs of time. 

As one aspect of the present invention, one or more of the unknown variables for 
several satellites are estimated for at least two base stations, and preferably for at least 
three base stations by a process that is performed at the rover. Using the numbers 1, 2, 
and 3 to represent the base stations Bl, B2, and B3, respectively, in the difference 
operator " A qr " the estimated values of the between-base station unknowns for each 
satellite 'V at each epoch are: 



Base Stations 1 and 2 - 


A 2,l*k > 




A 2A N L1 > S , 


A 2>l N L2 ' s 


[6A] 


Base Stations 1 and 3 - 


A 3,l T k> 




A 3>1 N L1 ' S , 


A 3A N L2 >* 


[6B] 


Base Stations 2 and 3 - 


A 3,2*jfc > 


H21I 


A 3 , 2 N L1 >\ 


A 3y2 N L2 > s 


[6C] 



An exemplary estimation process is outlined in a section below entitled " Between Base 
Station Processing, Part II. " As one option, the estimation process may break down 
the fixed integer ambiguities into their components as follows: 
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Base Stations 1 and 2 - A 2X N LX ' S , A 2l i// L1 (from A 2l N L1 ' s ) 

A 2A N L2 > S , A 2AW L2 (from A 2>l N L2 ' s ) 

Base Stations 1 and 3 - A 3l N L1 ' s , A 3l y/ LX (from A 3l N L1 ' s ), 

A XI N L2 ' S , A 3i# " (from A 3A N L2 > S ) 

Base Stations 2 and 3 - A 3 2 N L1 ' S , A 3 2 y/ LX (from A 3 2 N Ll,s ), 

H2 nL2 ' S > A 3)2 ^ L2 (from A 3 2 N L2,S ) 
In some applications of the present invention where lower estimate accuracy is 
acceptable, floating-point ambiguities maybe used in place of the fixed-integer 
ambiguities, and the following set of between-base station unknowns may be generated: 



Base Stations 1 and 2 - 


" A 2,l^k> 


A 2,l'*> 


A 2>l N Ll >\ 


A 2A N L2 ' S 


[7A] 


Base Stations 1 and 3 - 


- A 3,l*"jfc, 


A 3,l'*> 


A 3A N L1 '\ 


A 3>l N L2 >* 


[7B] 


Base Stations 2 and 3 - 


" A 3,2 T k> 


A 3,2^, 


A 3a N Ll >\ 


A 3a N L2 > s 


[7C] 



Between Rover and Base Station Processing. 

The processing of differences between the rover station and any one of the base 
stations is more complex because we do not know the position of the rover beforehand. 

Thus, the term A^ r i?£ = R* k -R s qk is not exactly known. However, given an initial 

estimated value of the rover's location, which we will denote as X 0 and which can be 
generated from the 50 Hz information signals by means known to the art, we can 
generate an initial range estimate between satellite "s n and the rover, which we denote as 

^0,* 9 311(1 we can g ener ate a linearized approximation of the term as follows, where we 
use r=0 to denote the rover and use q to generally indicate any one of the base stations: 

A *,o^ = K,k - = fe + 4 ■ ac\- R s q , k 
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where SX k is the vector difference between the actual location^ of the Rover and the 
initial estimated location X 0 (X 0 = X 0 + SX k ), where a s k is a row vector comprised of 
the partial derivatives of the distance form R s Qk =\\ X s - Xo || with respect to the rover 
location Xo, as evaluated at X 0 , and where A^ 0 #J k has been used as a shorthand 
notation for (ffj^ - R* k ). Formally, row vector a| is expressed as: 

The values of a s k depend upon the satellite location, which changes with each time 

moment and upon the initial location X 0 of the Rover, which may be periodically 

updated, particularly if the Rover is moving. The row vectors a s k for all of the satellites 

may be collected together as a single s-by-3 matrix, which we denote as A k . Matrix A k is 
commonly known as the Jacobian matrix, the geometry matrix, and the matrix of 
directional cosines. It is commonly computed in the art and described in many tutorial 
textbooks on global positioning. The reader unfamiliar with the GPS is directed to these 
tutorial texts for more detailed explanations of how matrix A k (which is oftentimes 
referred to matrix H or matrix G by these texts) is generated. 

With & q fiR s k = &qfi R Q y k + a l ' 8X , we may write the following versions of 

forms [4] (residuals) for the baselines between the Rover ("0") and any one of the base 
stations ("q yr ): 

*<lfi b i U = *qJ>Pk U ~ A <7 A" " AqfiTl [8A] 

A <7,0*f '* - A^pf 2 - W** - A *0** [8B] 

* qfi pi Us = a,, 0 ^ m - -±j A g , 0 R s k - A f >0 r; [8C] 

and the following versions of forms [5] 
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A *0*A = "i ' " cA <7,0** + A^o/f + A,,o^ M [9A] 



* q ,obi 2 > s =a s k -SX k -cA q90 T k + 



[9C] 



(/l zl >*j [9D] 

The solvable unknowns for form [9] includes SX k , as well as A qQ r k , A q0 I s k , 

&qfiN Ll > s , and A q0 N L2 > s . As described below in greater detail, estimation processes 

for these unknowns generally include instances of forms [9A]-[9D] for several satellites, 
and for several epochs of time. 

While the presentation of this section has been made with respect to a specific 
rover station (0), it may be appreciated that the discussion and above forms generally 
apply to any pair of receivers (q, r), where "r" substitutes for "0" in the above. 

Development of the Present Invention 

As part of making their invention, the inventors have recognized that the base 
station time-offset unknowns, when accurately estimated, should satisfy the following 
relationship: 

A 2 ,l^ + A 3a r k = A 3A r k . [10A] 
An underlying meaning of this relationship can be seen by noting that A 3 x r k = -A x 3 r k , 
and by rewriting the relationship in the equivalent form of A l x r k + A 31 r k + A x 3 r k = 
0. That is, the sum of the time offsets around a loop of base stations is equal to zero. 
The other unknowns, when accurately estimated, should satisfy the following similar 
relationships: 

*2j/f + A 3,2^ = A 3 ,,/| [10B] 
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A 2A N Ll > s + A 3>2 N L1 ' S = A xl N Ll ' s [IOC] 

A 2A N L2 ' S + A 3>2 N L2 >* = A 3A N L2 > S [lOD] 

Ai,W LX + A 3>2 y/ Ll = A^y/ Ll [10E] 

A u y/ L1 + A X2V L2 = A 3>x y/ L1 [10F] 

A 2>1 A>^ + A 3>2 N L1 >* = A 3A N L1 > S [10G] 

A 2>1 A>">* + A 3>2 N L2 ' S = A^ 2 -* [lOH] 

A 2)1 jV Z1 >* + A 3>2 N Ll > s = A^jV 11 ' 5 [101] 

A 2>l N L2 > s + A X2 N L1 ' S = A 3yl N L2 > s [10J] 



As a further part of making their invention, the inventors have recognized that the above 
relationships [10A]-[10J] should be satisfied by any set of three receivers, including sets 
where one of the receivers is the receiver of the rover station. We illustrate this point by 
replacing station B3 (3) with the rover station, which we will identify by the subscript 
number "0": 



A 2,in + A 0}2 r k 




[11 A] 


A 2 ,l/| + A 0>2 /* 


- A o,i/| 


[11B] 


A 2A N Ll > s + A 0 , 2 N Ll > s 


= A 0>l N Li > s 


[HC] 


A 2tl N L2 >° + A 0>2 N L2 >* 


= A 0A N L2 ' S 


[11D] 


A 2>l y/ LX + A^ 2 y/ Ll 


= *o,iv L1 


[HE] 


A 2 ,\¥ L1 +A 0)2 ^ Z2 


= A 0jX¥ L2 


[11F] 


A 2A N Ll > s + A 0i2 N Ll > s 


= A 0 ,iA> L1 ^ 


[HG] 


A 2>l N L2 >* + A 0>2 N L2 >* 


= Ao,!^ 2 '* 


[11H] 
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A 2 , X N LX > S + A 0 , 2 A^ = A 0 ,lN Ll ' S [HI] 
A 2A N L2 > S + A 0a N L2 > s = Ao 9 iN L2 > s [11 J] 

Still further, the inventors have recognized that relationships [10] and [1 1] can be 
generalized to apply to any loop of four or more receivers. Nonetheless, we provide 
exemplary embodiments of the present invention using loops of three receivers in order 
to simplify the presentation of the invention. 

According to general aspects of the present invention, the position of the Rover is 
estimated by forming a primary baseline between the Rover and one of the base stations, 
usually the closest base station, and then forming one or more secondary baselines from 
the Rover to other base stations. Then, one or more of the above relationships are 
applied to the configuration of stations to relate the measured data associated with the 
secondary baseline(s) with the measured data associated primary baseline using one or 
more of the unknowns associated with the baselines between base stations through forms 
[11]. To simplify the presentation of the present invention, and without loss of 
generality, we will use base station Bl to form the primary baseline with the Rover, and 
we use base stations B2 and B3 to form the secondary baselines with the Rover. Varying 
amounts measured data from the secondary baselines may be related with the primary 
baseline. In general, the greater amount of information so related increases the accuracy 
of the estimated position of the rover, and/or enables a greater spacing between base 
stations. 

First General Group of Embodiments 

In a first exemplary embodiment, the measured pseudo-range data from the 
secondary baselines are related with the measured pseudo-range data of the primary 
baseline. At first, only the LI -band data is used (forms [4A], [5 A], [8A], and [9A]). We 
take the case where the distances between the stations are small enough that the 

difference terms associated with the ionosphere terms ( A^ r /£ ) are negligible. For the 

primary baseline, forms [8 A] and [9 A] reduce to: 

A lfibt l ' S = Ai, 0 ^ U " A h0 R S k - A lj0 r; [12A] 
A lfibi hS =4 SX k -<*1,0** +A 1>0 ^ U [13A] 
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For the secondary baseline to the second base station, we have: 



*2,0bk U = A 2,op£ hS " A 2f0 £/ - A 2)0 7^ [14A] 
A 2i0 b^ s =a s k SX k -cA 2fi r k +A 2> o»f [15 A] 

For the secondary baseline to the third base station, we have: 

A 3 ,0^ U = A 3 ,oPf '* - A 3 , 0 ^ - A 3j0 ^ [16A] 

A 3,0b£ US =4 SX k - cA 3 9 0H +A 39 o"f ' 5 [17A] 



The terms in forms [12A], [14A], and [16A] are the residuals of differential navigation 
equations and are known. The unknowns are contained in right-hand sides of forms 
[13A], [15A], and [17A]. There are a total of six solvable unknowns: A xo r k , A 2 flT k > 

A 3,0 T k > the thr ^e components of SX k . However, using form [1 1 A], A 2 $T k ma Y be 

related to A 0 l r k through the estimated base-station time offset A 2 \T k as follows: 

A 2,0 T k = A 2j ir jt -A 0)1 r yt = A 21 r^ + A 1)0 r A: . In a similar manner, A 30 r k may be related 

to A 0 l z k through the estimated base-station time offset A 3l r k as follows: A 30 r k = 

A 3,\ T k- A 0 9 l T k = A 3jl z-£ + A lQ r k . These two forms reduce number of true solvable 

unknowns to four, and the following modified set of forms may be used to estimate the 
true solvable unknowns: 

A U<fik U =4 - 5X k -cA h0 T k +A l0 nj; l ' s [13A] 

(a 2) 0^ U +cA 2A r k )=a s k SX k -cA l>0 T k +A 2fi n^ s [15 A*] 
(a 3 ,0^ U +cA 3yl r k )=a s k -5X k -cA lfi T k + A 3j0 /if [17A*] 

As an equivalent, one may view forms A 20 T k = A 2) iT it + A lj0 r A . and A 3Q r k = 
&3,l T k + A\fiT k as increasing the total number of forms by two, and may then estimate 
the six unknowns the expanded form set: 

Al,0*f M = 4 ■ 5X k - cA h0 T k + A lfi n^ s [13A] 
A 2 ,0^ U =4 SX k -cA 2fi r k +A 2>0 n^ s [15A] 
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A 3) 0^ U = 4 ■ SX k - cA 3fi r k + A 3fi n^ s [17A] 
A 2 ,o^ = A 2( ir / t + A 1)0 ^ [11 A*] 

A 3,0T/c = A^r^ + Ajo^ [11 A"] 

The rover's location (X 0 = ^ 0 + ^ ) ma y be estir nated from each of the above 
form sets using any one of several know methods. Here, we demonstrate estimating the 
rover's location with a least squares fitting approach based on the first form set. After 

generating the residuals A Xfi b^ s for N satellites (s=l to s=N) for the baseline between 
the rover and the first base station (R-Bl), we collect them together into a first vector: 
A,,oBf =[A,, 0 Af.', AuBfU .... A^f. 

We do the same for the other two baselines, but modify their residuals according to the 
between-base station time offsets cA 2} \Tj c and cA 3 ^ : 



a n* Li - 
A 2> oB /t - 



[(A 2 ,0^ ' + cA 2,l^X (A^O^+cA^), .... (Ho^ hN +cA 2il T k )] T 

a 3 ,ob; zi = 

[(As.O^ + cA^MAs^ 

where the superscript " * " is used to indicate the modification. The noise terms can be 
similarly grouped as: A^n^ 1 , A 2>0 ii£ 1 ,and A 3 0 n^ 1 . The first form set of [13A], 
[15 A*] and [17A*] may then be written in matrix form as: 



A l,0Bf " 




A * 


-i 


a n* L1 




A * 


-i 


A R*^ 1 
A 3,0 B k 




A * 


-1 



A i,o«f 

A 2)0 n A 
A 3,onf 



[18] 



The solvable unknowns 
follows: 



cA l,0 T k 



may be estimated by a least squares process as 
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-1 



•(aaF-c- 1 *- 



cA l,0 T k 




[19] 



where: 




-1 



-1 



-1 



and where C n k is a 3N-by-JN covariance matrix for the noise terms A 10 n^ 1 , A 2 o^ 1 , 
and A 3 0 n^ 1 . The generation of covariance matrix C n fc is well known to the GPS art 



and described in the technical and patent literature (see for example U.S. patent No. 
6,268,824, incorporated herein by reference, where the covariance matrices are denoted 
as matrices R), and a description thereof is not needed for one of ordinary skill in the 
GPS art to make and use the present invention. (In a subsequent section, we describe an 
estimation process that we prefer to use, and which may be applied to these 
embodiments.) 

For the reader not familiar with the GPS art, we briefly note that matrix C n k 

generally comprises a diagonal matrix, with each diagonal element being related to the 
noise sources in the two receivers that define the baseline (in this case the rover and one 
of the base stations). A covariance factor is usually associated with each satellite signal 
received by each receiver, and this covariance factor is usually related to the signal-to- 
noise ratio of the signal (as received by the receiver) and the elevation angle of the 
satellite (the multipath error has a strong correlation with the elevation angle). Each 

diagonal entry of matrix C fc usually comprises an addition of the two covariance 
factors associated with the two receivers which contributed to the underlying noise 
quantity, e.g., base station Bl and the rover for noise quantity A l 0 n% l,s . For more 
details and information on the generation of the covariance matrices, the reader is direct 
to A. Leick, GPS Satellite Surveying, John Wiley & Sons. (1 995). Matrix A* is known 
as the observation matrix, and it relates the solvable unknowns to the residuals. The 
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solvable unknowns may also be estimated by other processes, such as various Kalman 
filtering processes. 

The above may be carried out by omitting the rows related to one of the 
secondary base stations B2 and B3, such as the last N rows associated with the third base 
station. The above may also be carried out by adding rows related to a fourth base 
station (and even more base stations). 

The second form set comprising [13A], [15A], [17A], [1 1 A*], and [1 1 A**] may 
be written as: 



Wi 1 
_ a _¥_b?_ 





I 


0 


0 


A *l 


0 


-1 


6 




0 


0 


-l 


o ! 


-1 


1 


0 


o i 


-1 


0 


1 



_sx k _ 

cA i,o T k 

cA 2,0 T k 



+ 



LI 



A «^1 
A n Ll 



0 
0 



[20] 



and the solvable unknowns may similarly be estimated by various least squares 
processes, and Kalman filtering processes. The observation matrix in the above form is 
denoted as A^* . 

Thus, the above exemplary embodiments provide methods of estimating the 
location of the rover station (R) with the use of a first base station (Bl) a second base 
station (B2), and optionally a third base station (B3) or more base stations. In summary, 
the locations of the base stations were obtained, and one of the base stations (e.g., Bl) 
was selected to form a primary baseline with the Rover station. We refer to this base 
station as the primary base station, and the other base stations as the secondary base 
stations. Additionally, for each of the secondary base stations, the time offset 
representative of the time difference between the clocks of the primary station and 
secondary base station is obtained. Also, measured satellite data as received by the 
rover, the primary base station, and the secondary base station(s) is obtained. From this, 

the set of residuals A^qB^ 1 associated with the primary baseline (R-Bl) is generated, 
and is related to the measured satellite data received by the rover station and the first 
base station, the locations of the satellites, and the locations of the rover station and the 

first base station. Similarly, the set(s) of residuals A 2) qB^ 1 , A 30 b£ 1 , etc. associated 
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with the secondary baseline(s) are generated, each set of residuals related to the 
measured satellite data received by the rover station and the secondary base station, the 
locations of the satellites, and the locations of the rover station and the secondary base 
station. Thereafter, the rover's location is estimated from the above sets of residuals, the 
time offset between the clocks of the base stations, and typically an observation matrix. 

The above embodiments can be equally implemented using the L2-band signals. 
The forms for the primary baseline (neglecting ionosphere terms) are: 



Ai,o^ 2 ' 5 = A i,oP^' S -A l>0 R s k - A 1>0 7jj 



L2,s 



A l,0 b k 2 ' S = a k 5X k ~ cA l,on + Ai,o"£ 



L2,s 



[12B] 
[13B] 



For the secondary baseline to the second base station, we have: 

A 2fi b£ 2 > s = A 2f0 p£ 2 > s - A 2fi R° - A 2fi T* [14B] 



A 2>0 b£ 2 > s =a s k SX k -cA 2>0 T k +A 2fi nj; 



L2,s 



[15B] 



For the secondary baseline to the third base station, we have: 

,L2,S — A _ - „ L2 ,S _A..B« A _ tS 



A^b^" =A X0P p>-A X0 R> k -A 



3fi T k 



[16B] 

Horf 2 '* =4 -SX k -cA 3fi T k +A X0 n^ s t 17B ] 
Not only may these forms be used in place of forms [12A]-[17B], but they may be used 
in addition to forms [12A]-[17B]. In the latter case, the first form set [18] expands to 



A, n B 



1,0% 
*L\ 

k 

*L\ 



A 2)0 B 



A 3,qB* 

A l,0Bf 

A Tt* L2 
A 2,0 B k 

A n* L2 



sx,. 



_ cA l,0 T k_ 



a 

A n Ll 
A 2,0"k 

A n L1 
A n i2 

a n L2 
a „£2 



[18 + ] 



where: 



A n*L2 _ 



K A 2,0b^ 2 ' l + cA 2A r k ), (A 2fi b^' 2 + cA 2>1 T k ), (A 2)0 ^ 2 ^ + cA 2)1 r /t )] 1 
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A 3,0« /t - 

[(Aa.O^' + cA^ ), (A 3>0 b£ 2 > 2 + cA 3>1 T k ),..., (A 3>0 ^ 2 '"+cA 3>1 r* )] T 

From form [18+], the solvable unknowns may be estimated by a least squares process. 

One might envision a different approach whereby the Rover's position is 
estimated by using each of the three baselines separately to generate three separate 
estimates of the Rover's location, and then averaging the three separate estimates to 
generate a final estimate. However, the present invention provides more accurate results 
than this possible approach because the present invention uses the additional information 
provided by forms [11] (specifically, [1 1 A], [1 1 A*] and [1 1 A**]), which correlates the 
values of the unknowns associated with the baselines. The achievement of more accurate 
results by the present invention is also true for the embodiments described below. 

Second General Group of Embodiments 

The second group of embodiments builds upon the first group of embodiments by 
generating and utilizing the residuals derived from the phase measurements at the 
receivers. We take the case where the distances between the stations are small enough 

that the difference terms associated with the ionosphere terms ( A^/^ ) are negligible. 

In later embodiments, we will take these ionosphere terms into account. With this 
assumption, forms [8C], [9C], [8D], and [9D] for the primary baseline become the 
following forms [12C], [13C], [12D], and [13D], respectively: 



W U = A lj0 ^ - A U0 Ri - A 1>0 7? [12C] 

Vt 1,1 =^h a l SX k ~f LU ^r k + Ay,** 1 -' +A li0 v^* [13C] 

A, 

Ai,o/>f = Ay,^* - -^j^i ~ -g&Wi [12D] 

A U>^ 2 ' 5 = ^l2j a l SX * -f L2 ' S *\,on +Al, 0 N L2 ' S +A 1>0 vf ■* [13D] 



For the secondary baseline to the second base station, we have: 
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A 2 ,0^ U = A 2j0 ^ U - ^ A 2i0 *f - ^ A 2>0 r; [14C] 

A W P F' S = jh a k dX * -f Ll ' SA 2,on +±2, 0 N L1 ' S +A 2>0 v^* [15C] 

A^jtf = A 2;0 ^ 2 ^ - A 2>0 ^ - -±- A2fiT s [14D] 

A 2,0Pf = ^7< -/ L2 ' 5 A 2)0 r, +A 2)0 ^"' S +A 2>0 vf •* [15D] 

For the secondary baseline to the third base station, we have: 

HOP?'* = A 3 ,otf'* " A 3 , 0 ^ " ^7 A 3)0 r/ [16C] 

A 3.0^ U = ^7«* ^ -f^Lz&k +*l,oN LU + A 3fi v^ [17C] 

A W*' = A 3)0 ^ 2 ' S - pL_ A 3i0 ^ - -±- A 3 , 0 ^ [16D] 

A 3.0^ 2 '' = T^j4 -/"''A^, + A 3)0 A^ 2 >* +A 3>0 v^* [17D] 

The terms in forms [12C, D], [14C, D], and [16C, D] are the residuals of the differential 
navigation equations and are known. The unknowns are contained on the right-hand 
sides of forms [13C, D], [15C, D], and [17C, D]. For Af satellites, there are 6*N 
equations and a total of (6+ 6*N) solvable unknowns: A l0 r k , A 20 T k > A 3,0 T k > the three 

components of SX k , and N instances of each of A X qN LI > s , A 2Q N Ll > s , A 30 N Ll > s , 

A l0 N L2 > s , A 2 flN L2 ' s , and A 3fi N L2 > s for s = 1 to s = N satellites. However, using 
forms [1 1 A], [1 1C], and [1 ID], the number of unknowns may be reduced to 4+2*W. As 
we saw above, A 2fi r k may be related to A 0 ^z k through the estimated base-station time 
offset A 2x r k using form [1 1 A] as follows: A 20 r k = A 2 jT k -A 0i iT k =A 2jl T k + A li0 T k . 
In a similar manner, A 30 r k was related to A 0 ir k through the estimated base-station 
time offset A 3A r k as follows: A 3fi r k = A 3)1 r^-A 0> i^ = A 3>1 r* + A Xfi r k . Form[llC] 
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fL\,s 
ktLU 



may be used to related each of the LI band ambiguities of the secondary baselines to the 
primary baseline as follows: 

A 3fi N Ll > s = A XI N LI > S - A 0>l N L1 > s = A 3A N Ll >° + A^N 1 . 

5 Similarly, Form [1 ID] may be used to relate each of the L2 band ambiguities of the 
secondary baselines to the primary baseline as follows: 

&2fiN L2 > s = A 2A N L2 > S - A 09l N L2 > s = A 2)1 A^ 2 '*+ A h0 N L2 > s 

A 3j0 ^ 2 ^= A 3jl ^ 2 ^ - A 0)1 ^ 2 ^ = A 3jl A^ + A hQ N L2 > s . 
Forms [15C] and [15D] for the second base station B2 may then be modified as follows: 

10 X ' [15C*] 

+ A 1>0 ^ + A 2>0 vf^ 

(a 2> 0^ +f L2 ' S A 2 , l r k -^N L ^)=-±-al .SX-f^A l0H 

X ' [15D*] 

+ A Xfi N L2 > s + A 2;0 vf'* 
Forms [17C] and [17D] for the third base station B3 may be similarly modified: 

A [17C*] 
■fA 1>0 ^ + A3, 0 v^ 

b,0P L k 2 > S +f L2 ' S A 3A r k -±3,l* L2 ' S )=-±ja S k -SX-f^A^n 

* r ' s . [17D*] 

+ A 1> o^^ + A 3> ovf^ 

15 Similar to form [18 + ], forms [13C, D] and modified forms [15C*, D*] and [17C*, D*] 
can be written in matrix form: 
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A 1,0P A 
A 2,0P^ 

A i,oP* 

A 2,0P* 
A n *L2 
A 3,0P* 



1 

AZ} 


•A, 


' T 1 

1 _ ^Ll 


! Jawv 


! °AfrAf 


1 

azI 


•A, 


1 T 1 
I — f 


! l NxN 


] °A^A^ 


A^l 


"A, 




| I A^7V 




AZ2 


•A* 


! _f^ 2 


! Q NxN 


! l NxN 


AZ2 


•A* 




! °AtaV 


! 'a^a^ 


AZz 


•A* 




°NxN 


' 1 NxN 



A,, 0 N 



L2 



V?- 

A 2,0^ 
A 3.0^ £1 

A i.o"* 2 

A 2 )0 vf 
A 3,0vf _ 



[21] 



where: 



• a 1,op£ 1 is an AT x 1 column vector of the LI -band phase residuals A l 0 p% l * s 
associated with the primary baseline, 

• A 2,0P*k Ll is an A^x 1 column vector of the modified LI -band phase residuals 

( A 2,0/>* ^ + f Ll ' s A 2i iT k - A 2y \N LlyS ) associated with the secondary baseline 
between the rover and base station B2, 

• A 3j0 p k is an Af x 1 column vector of the modified LI -band phase residuals 

[ A 3,0Pk l ' S + f Ll,S A 3,l r & ~ A 3,\N Ll ' s ) associated with the secondary baseline 
between the rover and base station B3, 

• A 1,0P* 2 is an A^x 1 column vector of the L2-band phase residuals A l0 p% 2 > s 

associated with the primary baseline, 
*L2 • 

• A 2) oP £ isanNxl column vector of the modified L2-band phase residuals 

{ A 2,0Pk 2,S + f L2,SA 2,l T k ~ &2 9 lN L2 ' S ) associated with the secondary baseline 
between the rover and base station B2, 
*L2 - 

• A 3j oP k is an N x 1 column vector of the modified L2-band phase residuals 

( A 3,0/>£ 2,5 + f L2 ' S A 3,l r * ~ A 3,lN Ll9S ) associated with the secondary baseline 
between the rover and base station B3, 

• is the derivative matrix previously described, 

• A Ll is a diagonal matrix of the carrier wavelengths of the satellites in the Ll- 
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band, and is its inverse matrix, 

A L2 is a diagonal matrix of the carrier wavelengths of the satellites in the L2- 
band, and A i\ is its inverse matrix, 

f Ll is a column vector of the carrier frequencies of the satellites in the LI -band, 

where A n - f Ll = c • 1 , where c is the speed of light and 1 is a column vector of 

l's, 
L2 

f is a column vector of the carrier frequencies of the satellites in the L2-band, 



where A 



v L2 * = 

• I js/xN is an TV x N identity matrix, and 

• 0 NxN is an N x TV matrix of zeros. 

The solvable unknowns \sX k , A 10 r^, A 10 N L1 , A^N^J^may be estimated by 
applying a number of processes to the combination of forms [21] and [18+]. Further 
below, in a subsequent section, we described a preferred process that can be used on the 
combination of forms [21] and [18+]. Here, we describe how a least squares process 
may be applied to combined forms of [18+] and [21]. 

To combine form [21] with [18+], we rewrite form [22] to substitute c • A l 0 r k 

for A 1)0 t-£ as one of the unknowns, using the facts that A n 



•i Ll =c l and 



L Z,2 



>f L2 =c-i: 



A i,oP* 

A 2,0P* 
A 

A 3,0P* 
A n L2 

a *L2 
A 2,0P* 

A ^*£2 
L A 3,0P* J 



All 


• A * 


l- A zl 


1 


! I NxN 


1 °NxN 


All 


" A * 


:-az! 


1 


l NxN 


j °NxN 


AZl 


• A * 


i-AZi 


1 


l NxN 


! ®NxN 


az!> 


" A * 


|-AZ^ 


1 


°NxN 


! * NxN 


AZ^ 


■ A k 


"AZ^ 


1 


®NxN 


! I NxN 


AZ^ 


■ A k 


-AZ2 


1 


°NxN 


! 1 NxN 



5X, 



?: A i.Pf* 

" L2 



Ai >0 N J 



+ 



A 2,0^ 1 

A 2,0'f 
L A 3,0vf . 



[22] 



This facilitates the combination of forms [18+] and [22] as: 
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A R*^ 1 

A 2,0 B A: 

A 3,0«A: 




A i,oP* 

A 2,0P k 

a n *L2 
L A 3,0P* J 





! -i 








K 


! -I 








A k 


-l 








A k 


| -l 








A k 


1 -i 




i ®3NxN 


! O^Wv-AA 


A k 


i -I 








A l\' A k 


:-az! 


i 






AZ}-A^ 


I-azI 


i 




! ®3NxN 


AZi'Aj 


!-az! 


I 






A Z2-A, 


"AZ^ 


• I 






A Z2'A A 


-AZ2 


i 


°3NxN 


l 3NxN 


AZ^-A, 


~ A L2 


I 







8X, 





[23] 



A least squares process is applied to form [23] to generate the floating ambiguities over 
several epochs, and the floating ambiguities are averaged to generate a estimate of the 
floating ambiguities. For each epoch, the following may be generated: 





sx k 






5 


°- A \,o T k 
Ai, 0 N L1 








_A l5 oN i2 







g^LX 



0 



Q 3Nx3N °3Nx3N °3Nx3N 
L2 

] 3Nx3N C n,k ®3Nx3N °3Nx3N 

°3Nx3N 



-i-l 



°3Nx3N °3Nx3N 
°3Nx3N ®3Nx3N °3Nx3N 



r*L2 



A 



23,k 



^n,k 


®3Nx3N 


®3Nx3N 


°3Nx3N 


3Nx3N 




®3Nx3N 


°3Nx3N 


3Nx3N 


°3Nx3N 




°3Nx3N 


3Nx3N 


°3Nx3N 


°3Nx3N 





-1-1 



A l,oBf 

A R*^ 1 

a n* L1 

A 3,0 B /fc 




A n L2 

a ~*L2 

a n *L2 
A 3,0P* 



-1 
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where A 2 j^ is the observation matrix of form [23] (the matrix which is multiplying the 
solvable unknowns in form [23]), and where different covariance matrixes Cy] k and 

C yk are used for the LI -band and L2-band data since the magnitudes and variances of 

the noise sources associated with the phase measurements are different from those of the 
noise sources associated with the pseudorange measurements. Also, the least squares 

method normally produces the floating ambiguities A 10 N Z1 , A^qN^ 2 rather than the 

fixed-integer ambiguities A^qN^ 1 , A 10 N l2 . In practice, the least squares process is 
applied over many epochs, and the computed floating ambiguities are averaged to 
generate a final estimate of the floating ambiguities. Several averaging processes are 
described in U.S. patent No. 6,268,824, which is incorporated herein by reference, and 
may be used. After generating a suitable set of floating ambiguities, a conventional 
method of generating the fixed-integer ambiguities or the integer ambiguities may take 
place. 

With the fixed-integer ambiguities estimated, the unknowns SX k and c • A 10 r^ 

may be estimated by substituting the estimated values of A 10 N Z1 and Aj qN 1,2 into 

form [23], and moving these terms to the left-hand side with the residuals to provide 
form [24]: 




} 3Nx\ 




L2 



A * 


-1 




A * 


i -I 




A * 


! -l 




A * 


-i 




A * 


-i 




A * 


i -i 






j-AZi 


1 


A Zi ,A * 


!-AZ| 


1 


AZl-A^ 


i-Ai! 


1 


A l2" A * 


!- AZ2 


1 


A l2" A * 


-AZ2 


1 


A £2 'A A 


-AZ2 


1 



sx, 



_ c - A l,O r yfc. 



A l,0 n * 
A 2,0» k 
A 3,0 n £ 




L2 



A 2,0^ 

a L2 
A 3,0^ 



[24] 
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We consider the left-hand side of form [24] to be a plurality of sets of residuals since 
each comprises the known quantities of differential navigation equations. A second least 
squares process may then be applied based on form [24] to estimate the unknowns 

8X k and c-A^qt^ as follows: 



SX, 





°3Nx3N 


®3Nx3N 


0 


3Nx3N 


^n 9 k 


°3Nx3N 


0 


3Nx3N 


°3Nx3N 
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3Nx3N 


°3Nx3N 


°3Nx3N 


( 



3Nx3N 
3Nx3N 
3Nx3N 
'v,k 



-1 



'24,k 



^n,k 


°3Nx3N 


°3Nx3N 


°3Nx3N 


3Nx3N 


^n,k 


°3Nx3N 


°3Nx3N 


3Nx3N 


®3Nx3N 


*~y,k 


°3Nx3N 


3Nx3N 


°3Nx3N 


°3Nx3N 


^v,k 



-1 
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a 3) ob;- 




°3Nxl 




A l,0Bf 








a n* L2 




°3Nxl 




A n* L2 










A l,oNf ' 

A 2fiK 
A 3, 0 Nf 




A n L\ 
A 1,0PA; 

A 2,0P,fc 
A 3,0P^ 








A l,0Nf 
A 2,0Nf 
A 3,0Nf J 




A i,oPf 

A 2,0P k 
A n *L2 





[25] 



3fC 

where A 2 4 9 k is the observation matrix of form [24] (the matrix which is multiplying the 

solvable unknowns in form [24]). In form [25], the floating ambiguities may be used in 
place of the fixed-integer ambiguities. However, lower accuracy generally results, 
although the estimation speed is increased since the step of generating the fix-integer 
ambiguities from the floating ambiguities may be omitted. 

With large amount of measurement data afforded by the present invention, 
various residuals may be omitted from use in the estimation process. For example, we 
may only work with LI -band data of all three stations and only use the residuals 
associated with this data. We may also only work with data and residuals from two base 
stations (the primary baseline and one secondary baseline). For applications needing 
lower accuracy, we may work with the phase and pseudorange data of the primary 
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baseline (R-Bl) and just the pseudorange data or phase data of one of the secondary 
baselines. Also, one can undertake an analysis of the satellite constellation and select the 
satellites which should provide the highest accuracy, using the pseudo-range and phase 
data from the LI - and L2- bands. Finally, while we have illustrated the invention with 
the single-difference navigation equations, it may be appreciated that higher-order 
differences of the navigation equations may be used. Such higher-order differences have 
known quantities (which form the residuals) and similar unknowns to be solved for, 
which can be solved by the above-described methods. 

Thus, the above exemplary embodiments from the second group of embodiments 
provide methods of estimating the location of the rover station (R) with the use of a first 
base station (Bl) a second base station (B2), and optionally a third base station (B3) or 
more base stations. As a summary of an exemplary embodiment of the second group, the 
locations of the base stations were obtained, and one of the base stations (e.g., Bl) was 
selected to form a primary baseline with the Rover station. Additionally, for each of the 
secondary base stations, the time offset representative of the time difference between the 
clocks of the primary and secondary base stations is obtained, and a set of satellite-phase 
cycle ambiguities related to the baseline between the secondary and primary base stations 
is obtained for both of the frequency bands. Also, measured satellite data as received by 
the rover, the primary base station, and the secondary base station(s) is obtained. From 

this, the sets of residuals A lj0 B^, A 10 B^ 2 , A 1>0 pf and A^p" of differential 
navigation equations associated with the primary baseline (R-Bl) 'are generated, and are 
related to the measured satellite data received by the rover station and the first base 
station, the locations of the satellites, and the locations of the rover station and the first 
base station. Similar set(s) of residuals associated with the secondary baseline(s) are 
generated, each set of residuals related to the measured satellite data received by the 
rover station and the secondary base station, the locations of the satellites, and the 
locations of the rover station and the secondary base station. Thereafter, the rover's 
location is estimated from the above sets of residuals, the time offset between the clocks 
of the base stations, the sets of satellite-phase cycle ambiguities related to the baseline 
between the primary base station and the secondary base stations, and typically an 
observation matrix. 
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Third Group of Embodiments 

In the above embodiments, the ionosphere delays were assumed to equally affect 
the rover and base stations, and the between-station differences were neglected. When 
the base stations and rover are separated by large distances, better accuracy may be 
5 obtained by taking into account the ionosphere delays. This can be done in a number of 
ways according to the present invention. FIG. 3 provides a representation of the 
ionosphere delays of one satellite "s" at the base and rover stations. Shown is a 3-d 
Cartesian system having two planar axes, north (n) and east (e), to represent the terrain 
on which the stations are located, and a vertical axis to represent the ionosphere delays of 
10 a satellite "s" as a function of the terrain. The ionosphere delay will be different for each 
satellite. The locations of the rover station R and three base stations Bl, B2, and B3 are 
indicated in the north-east plane of the figure. The ionosphere delays for each of these 

stations are indicated as I$ k , l{ k , and I\k > respectively. The preferred 
embodiments of the processing of the base station data, which is more fully described 
15 below, generates estimates for the between base stations differences in the ionosphere 

delays A 12 /| , A 13 /| , and A 2> 3/£ , which we call ionosphere delay differentials. From 

two of these differentials, an estimate for the ionosphere delay between the rover and any 
of the base stations can be generated. Here, we are most interested in the difference 

^Ifl^k associated with the primary baseline, the one between the rover and the first base 
20 station. 

The first approach is to generate an estimate A 1>0 /£ °f ^\Mk based on an 
interpolation of two of the known ionosphere delay differentials ( A 1>2 /| , Aj 3 /| , 
A 2,3^l ) onto the known approximate location of the rover. Here, we use A 12 /| and 
Ai^/jt in the following interpolation: 

25 

Ai j0 7/=«Ai, 2 /|+>a\ u /| [26] 

where a and /? are interpolation constants. The interpolation constants are determined as 
follows. Let Xj 9 X2, and X3 be three-dimensional vectors which represent the positions of 
30 base stations Bl, B2, and B3, respectively. Let Xq^ represent the approximate 
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estimated position of the rover at the "k" epoch. Furthermore, let the notation {X} n 
denote the north component of a position vector X or a difference vector X of positions, 
and the notation {X} e denote the east component of a position vector X or a difference 
vector of positions X. We can then write the north and east components of the position 
5 differences {X ok -X l } n and {X 0k -X x } e as: 



There are two equations in the two unknowns a and /3, which can be readily solved for a 
and 0. When Xq^ is within the triangle formed by the base stations, both of a and /3 are 
10 greater than zero, and their sum is equal or less than one. 



ionosphere delay does not always vary in a linear manner over the terrain, and often has 
second order variations with respect to the east and north directions, which are not 



20 which is graphically shown in FIG. 3. 

With form [28], estimates for ionosphere differentials on the secondary baselines 

A 2> o^ and A 3)0 /£ can be generated using form [1 IB] (A 2> i/£ + A 0,2 7 * = A 0,l 7 ^) as 
follows: 



{Xo,k ~X l } e =a{X 2 -X x } e + /3{X 3 -X x } e 



n 



[27] 



Because X o k is not necessarily at the exact position of the rover, the estimate 
A l,0^fc S wi U not necessarily be equal to the true value A 1)0 /£ . In addition, the 




A 2) 0'jfc = A 2)1 /| + Ay,/* + S h0 I s k 



[29] 



25 



Wk = Wk + a 1>0 7/ + s l0 i s k 



[30] 
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In composing forms [12] - [17], we neglected the ionosphere delays. In the third 
group of embodiments, we add the ionosphere delay terms to these forms. By comparing 
the general forms [9A-D] to the specific forms [13A-D], [15A-D], and [17A-D], we can 
generate augmented forms [12]-[17] as follows: 
The primary baseline - 

Al.obJ' 1 = A^opI 1 -Aj.oR* -*i,oT k -AiA [12AT 

Ai, 0 b^ = A k -dX k -cA Xfi r k 1 + 8 lj0 l£ + A^n^ 1 [13AT 



Ai, 0 b* 2 = a i,op* " a i,oR^ -Ai, 0 t a -a1 2 All *i,oh [12BT 

Ai, 0 b" = A^ ■6X k -cA 1>0 r k -i + Al 2 All °1,<)I* + *1,0"¥ [13B1 

Ai, 0 P^ = A 1>0 <p" " A Ll A l,oR^ " A£i A^ +A£\ A uo l k [12C] 

A,,oP^ = A£\ -A k -8X k -f Ll -A h0 r k 1 

-A£i 8 1 , 0 I jt+ A 1 , 0 N L1 +A 1>0 vJ' 1 



[13CT 



A 1,op]T = A l,0^ 2 - A L ! 2 A^qR^ - A L2 A h0 T k 

+ A L2 A L1 A l,0^t 
Ai,oP^ = Al} 2 -A k .8X k -f L2 .A 1>0 r k 1 +A 1)0 N L2 

-A L2 A L1 8 10 I^ +A 1>0 v jt 



[12D-] 



[13D"] 



The secondary baseline to the second base station B2- 

r,. -1 

[14A] 

Ai^b^ 1 = A* -dX k -cA l0 r k 1+ 5 lj0 I k + A 2>0 nJ' 1 [15A1 



*L1 LI rw 

A2,0 b /t = A 2 ,0PA: - A 2,0 R yfc ~ A 2,0 T )t + cA 2 j? k • I 

-(a2,i^+a 15 oTJ 



A 2,ob^ L2 = A 2 , 0 pJ' 2 -A 2 , 0 R A -A 2 ,o T A +cA 2 ,i^ I 



- Al 2 A l J (A 2fl l k +A h0 l k ) 



[14BT 
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A 2,0^ L2 = A A -8X k -cA l0 T k 1 

2 2 L2 t 15B ^ 

+ A£ 2 A L f Sj^+Az^nL 2 



^oP* = A^ 1 - A£i • A 2j0 R A - A£\ - A w T t 

ri4cn 

+ /" • A 2 ,ir A • 1 - A 2>1 N L1 + A£\ ■ (a 2)1 I* + A lj0 T t ) 

A 2 , 0 P^ L1 = A£\ -SX^ .A lj0 r A -1 

+ A 1>0 N L1 - A£j 8 1)0 I* + A 2)0 vL 1 



A 2 , 0 P* = A 2 , 0 <pP - A L2 • A 2 qR^ - A L2 • A 2 0 T, 

[14D1 

+ /" • A 2> i^ • I - A 24 N L2 + A L2 - A£? • (a 2j1 I* + A l5 oI* ) 
A 2 , 0 P^ L2 = A£ 2 • A fc - dX^ - f L1 ■ A h0 r k ■ 1 

[15DT 

+ A 1>0 N L2 -A L2 A£? 8 1)0 I^ +A 2>0 vP 



The secondary baseline to the third base station B3- 

*L1 LI — 

A 3,0*> k = A 3 ,0P^ " A 3,0 R * - ^3,0^ + cA 3A T k ■ 1 



-(a 3 ,iI a +a 15 oTJ 



[16A1 



A 3 , 0 b^ L1 = A k bX k -cA h0 T k 1 + 8 h0 l k +A 3>0 nJ' 1 [17A1 



A 3 ,obl L2 = A 3 ,oP^ 2 -A 3 , 0 R A ~^3,0 T k +cA 3fl r k 1 
~ a l 2 • A£ 2 -(a^Ij. +A l5 oT /t ) 

A 3 , 0 b^ L2 = A^ -8X k -cA l>0 r k 1 

+ Al 2 • A£ 2 -8 ly0 l k +A 3j0 n^ 2 



[16B1 



[17B1 



A 3,0P* 



A 3 , 0 <p" - A L \ A 



3f0 R k - A L \ • A^qT^ 



+ f Ll ■ A3.1T* • 1 - A 3>1 N L1 + A£\ • (a^I* + A 1)0 1* ) 



[16C1 
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-1 



Z,l 



k i,o 



N L1 -A£ 1 1 .5 l5 oI^A 3j0 v^ 1 



[17C1 



+ /" * ■ 1 - A 34 N L2 + A L2 • All ' faa 1 * + A lfih ) 

A 3,0P^ L2 = A L2 " A k ' - f L1 . Aj 0 rjt • 1 

+ Ai, 0 N L2 - A L2 • A£* .5^ + A 3s0 vP 



This set of forms comprises N additional solvable unknowns, as present in the vector 
§lfl l k > which have been incorporated into forms [13^ [15^ and [17]. Also, the rough 

ionosphere approximation A\$l k has been incorporated into forms [\2% [14], and [161. 

Finally, the between-base-station data A^il^ has been incorporated into forms [14*], and 

10 the between-base-station data A 3?1 I^ has been incorporated into forms [16"]. 

The above-described processes of estimating the floating ambiguities and 
thereafter estimating the rover position can be expanded to include the unknowns 
8 1,0 I A: • For lower accuracy, the unknowns b x ^l k may be omitted from the above 

forms [13], [IS*], and [17 1 ]. However, approximate information on the ionosphere 
1 5 delays is incorporated into forms [ 1 2 *] , [ 1 4 ] , and [ 1 6 *] . 

Embodiment Implementations 

Having thus described three general groups of embodiments, we now describe an 
exemplary rover station 100 in FIG. 5 that may be used to implement any of the above- 

20 described embodiments. Description of rover station 100 is provided in conjunction with 
a flow diagram shown in FIG. 6. Referring to FIG. 5, rover 100 comprises a GPS 
antenna 101 for receiving navigation satellite signals, and RF antenna 102 for receiving 
information from the base stations, a main processor 1 10, an instruction memory 1 12 and 
data memory 1 14 for processor 1 10, and a keyboard/display 1 15 for interfacing with a 

25 human user. Memories 1 12 and 1 14 may be separate, or difference sections of the same 
memory bank. Rover 100 further comprises a satellite-signal demodulator 120 for 
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generating the navigation data pfy s , pffi , #fy s , and <f>ffi for each epoch k from 

the signals received by GPS antenna 101, which is provided to processor 110. Rover 100 
also comprises a base-station information demodulator 130 that receives information 
signals from the base stations by way of RF antenna 102. Demodulators 120 and 130 
may be of any conventional design. The information received by demodulator 130 
includes the positions (Xi, X 2 , X 3 ) of the base stations, the satellite-navigational data 

(e.g., k, pffi , p^ s , <f>^ s , and fffi , r =1,2,3) received by each base station at each 
epoch k, and the between-base station unknowns (e.g., A r q r k , A r q N LU , A r y/ LX , etc. , 
{ r ><j} = {1> 2 }>{ 2 >3},{3,1}). Each set of information may be transmitted on arespective 
frequency channel. The between-base station unknowns may be generated by the first 
base station Bl, and thereafter transmitted to the rover by the first base station. The first 
base station may receive the satellite-navigation data from the other base stations in order 
to compute the between-base station unknowns. Methods of generating the between- 
base station unknowns are described below in greater detail in the section entitled 
Between Base Station Processing, Part II . As another approach, the between the 
station unknowns may be generated by rover 100 locally by a between-base station 
processor 140, which receives the positions of the base stations and their satellite 
navigation data from base station information to modulator 130. Processor 140 may 
implement the same methods described in greater detail in the below section entitled 
Between Base Station Processing. Part II . In addition, processor 140 may comprise its 
own instruction and data memory, or may be implemented as part of main processor 110, 
such as by being implemented as a sub-process executed by main processor 110. 

Main processor 110 may be configured to implement any and above described 
embodiments by the instructions stored in instruction memory 112. We describe 
implementation of these embodiments with respect to FIG. 6, where certain of the steps 
may be omitted when not needed by a particular embodiment. In step 202, the locations 
Xi, X 2 , X 3 of the base stations Bl, B2, B3 are received by base-station information 
modulator 130 and conveyed to main processor 110. These locations, and the location of 
the rover station, are measured at the phase centers of the GPS antennas. Thus, RF 
antenna 102 and demodulator 130 provides means for receiving the locations of the first 
base station and the second base station. Also in step 202, main processor 1 10 generates 
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an initial estimated location for rover 100, which is relatively course. This may be the 
center of the triangle formed by the base stations, or may be derived from a conventional 
single point GPS measurements (as opposed to a differential GPS measurement), or 
maybe generated by other means. The means for performing this step is provided by 
main processor 110 under the direction of an instruction set stored in memory 1 12. 

In step 204, one of the base stations is selected as a primary base station (Bl) to 
form a primary baseline with the Rover station. The selection may be arbitrary, or may 
be based upon which base station is the closest to the initial estimated location of the 
Rover. The other base stations are secondary base stations and form secondary baselines 
with the rover. The means for performing this step may be provided by the human user, 
as prompted by main processor 100 through keypad/display 1 15, or may be provided 
directly by main processor 110 under the direction of an instruction set stored in memory 
112. 

In step 206, main processor 1 10 is provided with the measured satellite- 
navigation data (e.g., p^ s , p ^ s , <f>Q s , and tffs ^ ag received by the rover 

the primary base station(r=l), and the secondary base station(s) (r=2,3) at one or more 
time moments k (epochs). The data is provided by the modulators 120 (for the Rover 
data) and 130 (for the base station data). Although these sets of data for each epoch k 
may be received at slightly different times (because the base stations are at different 
distances from the Rover), the data sets are time-stamped with the epoch identifier 
(which is conventional practice), and can be stored in a synchronized queue until all the 
data sets for the epoch are received. Also during this step, main processor 110 
determines the positions of the satellites for these time moments from orbital predictions, 
and generate the computed ranges of the satellites to the rover and base stations based on 
the positions of the satellites and the stations. The means for performing this step is 
provided by main processor 1 10 under the direction of an instruction set stored in 
memory 1 12, with the computed information being stored in data memory 114. 

In step 208, for each secondary base station, Rover 100 obtains a time offset 
representative of the time difference between the clocks of the primary and secondary 
base station (e.g., &2,\ r k , A 31 r^) at the time moments k. Rover 100 may obtain this 
information by directly receiving it from the primary base station Bl by way of 
demodulator 130, or may obtain at this information by using between-base station 
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processor 140 to generate it, as indicated above. Either of these approaches provides the 
means for obtaining these time differences. As used herein, the term "obtain" 
encompasses both the receiving of the information from an outside source (e.g., the 
primary base station) and the generating of the information locally by processor 140. 

The next step 210 is optional, depending upon the embodiment being 
implemented. In this step, rover 100 obtains, for each secondary base station, a set of 
satellite-phase cycle ambiguities related to the baseline between the secondary and 
primary base stations for one or more of the frequency bands (e.g., LI and L2) at the 

time moments k. These ambiguities may be in floating form (e.g, A 2 ^N Ll ' s > 

A 2yl N L2 > s , A^N 11 ' 5 , A 3 l N L2 > s ) 9 fixed-integer form (e.g, A^N 11 ' 5 , A 2fl N L2 > s , 

A 3,lN Ll,S > &\\N Ll9S ) or integer plus fractional phase form (e.g., A 2yl N LU , 

*VLV L1 l A 29l N L2 > s , A 2 ,iW L2 ; * 3 ,iN Ll > s 9 A 3ylV , Ll ; A xl N L2 >\ A 3 ^ L2 ). Rover 
100 may obtain this information by directly receiving it from the primary base station Bl 
by way of demodulator 130, or may obtain at this information by using between-base 
station processor 140 to generate it, as indicated above. Either of these approaches 
provides the means for obtaining this information. As used herein, the term "obtain" 
encompasses both the receiving of the information from an outside source (e.g., the 
primary base station) and the generating the information locally with processor 140. 

In step 212, which is optional depending upon the embodiment, main processor 

110 generates tropospheric difference terms A^ r 7£ at the time moments £, which are to 

be used in the residuals. In the step, main processor 110 also generates the observation 
matrices for each time moment k. The means for performing this step is provided by 
main processor 110 under the direction of an instruction set stored in memory 112, with 
the difference terms being stored in data memory 114. 

In step 214, main processor 110 generates the sets of residuals of the single- 
difference navigation equations (e.g., Aj qB^AjoB^ 2 , A 1)0 p^ and A 10 p£ 2 ) 

associated with the primary baseline (R-Bl) at the time moments k, and generates the 
sets of residuals of the single-difference navigation equations associated with the 
secondary baseline(s) (R-B2, R-B3) at the time moments k. The forms for generating 
these residuals were described above, and depend upon the embodiment being 
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implemented. The means for generating the residuals is provided by main processor 110 
under the direction of an instruction set stored in memory 112, with the residuals being 
stored in data memory 114. 

As an optional step 216, main processor 110 generates ionosphere corrections 
5 and adds corrections to the residuals. The means for performing this step is provided by 
main processor 110 under the direction of an instruction set stored in memory 1 12, with 
the corrections being stored in data memory 114. 

Finally, in step 218, main processor 110 estimates estimate the rover's location at 
the one or more time moments k from the sets of residuals, the time offset between the 
10 clocks of the base stations, the sets of satellite-phase cycle ambiguities related to the 

baseline between the primary base station and the secondary base stations (optional), and 
an observation matrix. The estimation may be done by the previously-described 
methods. The means for performing this estimation step is provided by main processor 
110 under the direction of an instruction set stored in memory 112, with the computed 
1 5 information being stored in data memory 114. 

Keypad/display 115 may be used to receive an instruction from the human user to 
commence an estimation of the position of the Rover, and to provide an indication of the 
estimated position of the Rover to the user. For some applications, it may be appreciated 
that human and interaction is not required and that the keyboard/display would be 
20 replaced by another interface component, as needed by the application. 

Preferred Floating Ambiguity Resolution Process 

We demonstrate the preferred floating ambiguity resolution process using forms 
[12] -[17'] above. It may be appreciated that the other previously-described 
25 embodiments of the present invention may use also this floating ambiguity resolution 
process by simply omitting forms and/or terms of forms from the process. The general 
view of the preferred ambiguity resolution process is to reduce the value of the following 
form during a series of epochs: 

30 =|(A 1 , 0 N it . 1 -A M N^) r D^ 1 (A ls0 N^_ 1 -A 1?0 N*) 
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+ ( A * 8X * - cA lOH l+8l,0l* - AqX^Ffe] -1 * ^ 



HZ 



3 [ ( A ,8X, -cA^r, -1+ A 2 2 A£j b lfi l k - A q ,^J [c^ 2 ]" 
[a^X* -cAj^ -1 + A^ A"5 S^I* -A q , 0 b* L2 ) 



1 ^ 

x 



q=\ 



q=\ 



( A L1 A * 8X * -/ L1 A 1>0 r A -T-A-l .a li0 I 4 + A l50 N LI -A^oP^flc^]" 1 x 
, ( A L1 A * 8X * -f L %,on * * - a li -8l,0l* +A 1)0 N L1 -A q , oP ; L1 ) 

( A L2 A * 8X A -/"Al,o^ l-A L2 A£2 6l0 i /t + A 1>0 N L2 - A qfi p^ 2 f[ C ^ ^ 



A A L2 A k* X k -f LlA lfi*k - l-A^A-^oI* +A 1>0 N L2 -A q , oP ; L2 
+ ^Sl,oI/t r [«l,o]" 2 8i, 0 I)t [31] 

Inform F £ , A^qN^ represent the combined vector of current estimates of floating 
ambiguities: A^qN* = [ A 1?0 N L1 , A 10 N L2 ], and the following are used for the cases of 
q=l: A q ,ob; L1 =A q>0 b^, A^ 2 =A q , 0 b^ 2 , A^p*" = A q , 0 pL 1 ,and 
A q ,oP^ L2 = A q5 oP^ 2 • We have also broken down the weighting matrices 
[ C n 1 Y ' [ c v 1 » eta. as follows: 



teiY 

0 


! o 


! o 

•1 T Q 


0 


\ ° 





0 


j 

> 0 


! o 
: ' 0 


0 


! o 





6tfc. All of the other terms used above have been previously described with the exception 
of covariance matrix a 1Q and matrix D. The between station ionosphere delay, 

generally denoted here as A^ r /£ , usually follows the Gauss-Markov time model: 
V 7 *+l = A *,r**> [GM1] 
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where y = exp) — - 



<"7 



, r e is the time between successive epochs [sec], Xj is the 

ionosphere model correlation time [sec], ts. q r s^ is white noise with the zero mean value 
and the variance of: 

[a£ r ] 2 =(l-r/)[o r £ r ] 2 > where [GM2] 

5 l^lrf =[a] 2 \\\x r -X q \ ] 2 10~ 12 [in dimension of m 2 ], 

where, or e [1,5] is a scaling parameter depending on the solar activity. Covariance 
matrix a 1Q is a diagonal matrix of all the individual values ct^q , and where the inverse 
square [<*i 5 o] of the covariance matrix is a diagonal matrix of the inverse squares 

[c^q] 2 of the individual values 5j 5 q. Since this process is only accounting for the 

10 second order effects in the ionosphere delay differentials (rather than the full amount), 
the value of a is generally 2 to 3 times smaller than the value used in a single baseline 
ambiguity resolution process which accounts for the full amount. 

Inform [31], the term ^(A^N^ - A^N* ) T D k -l(^lfi^k-l - A lj0 N*) is a 
cost function which effectively averages the floating ambiguities over several epochs by 
1 5 introducing a penalty if the set of estimated floating ambiguities A 1?0 N^ at the k-th 

epoch is too far different from the previous estimated set. Aj . With each 

successive epoch, the weighting matrix Dk-i is generally made to be more convex. The 
second term of F 2 generates the weighted sum of the squared residuals of forms [13A], 
[15 A], and [17A]. In a similar manner, the third term of F% generates the weighted 

20 sum of the squared residuals of forms [ 1 3B *\ , [ 1 5B ] , and [ 1 7B 1 , the fourth term 

generates the weighted sum of the squared residuals of forms [13C], [15C], and [17C], 
and the fifth term generates the weighted sum of the squared residuals of forms [130], 
[15D], and [17D]. The residual of each of the above forms of [13^, [15], and [17] is 
the difference between the right-hand side and left-hand side of the form. The weights 

25 are defined by the corresponding inverse covariance matrices. 

An exemplary estimation process for the floating ambiguities according to the 
present invention employs form [31] in an iterative manner. We start at the initial epoch 
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k=0 with the weighting matrix D 0 set to the zero matrix, and an initial guess of floating 
ambiguities A 1>0 N 0 equal to zero. The first term of form [31] evaluates to zero for this 
initial epoch. We then generate a set of values for 8X^ , A 10 r^ , 8 10 l£ , A^N* at a 
first epoch k=l which moves the value of F £ towards zero. This generates an initial 
estimate A^qI^ for the floating ambiguities and the rover's position (by way of 8Xj). 
A weighting matrix Di for the ambiguities is then generated, and a new set of values is 
generated for 8X* , A 10 r^ , 8 10 l£ , A 1?0 N^ at the next epoch k=2 which moves value 

of F s (k=2) towards zero. As a result, subsequent estimates A lj0 N 2 and 8X 2 are 
generated. This iteration process continues, with A 1>0 N^ and 8X k generally improving 
in accuracy as the iterations progress. 

Generating the values of 5X k , A lf0 r k , S^qI^ , A lj0 N£ to minimize form [31] 
can be done in the following manner. Form [31] is constructed in a manner whereby 
each term is generally of the form: (M-Y -b) T -W • (M-Y -b), where Y is a vector of 
unknowns (e.g., some or all of 6X k , A l 0 r k , 8 1>0 I^ , A^qN^ ), M is a matrix of 
constants which multiply the unknowns (e.g., A k , c), b is a vector of known values (e.g., 

A q,O b I L1 > A q,0 b £ L2 )> and W is a weighting matrix (e.g., D k , [c^ 1 ). The term (M-Y 

T 

-b) 1 W-(MY-b) can be reduced in value towards zero by generating a vector Y 
which satisfies the relationship (M T -W-M)-Y = M T -W-b, which we can simplify as 
H-Y = B. If H is not singular, we may apply an LU-decomposition method on H to find 
Y. The contributions of the six terms of form [31] may be synthesized according to the 
following form: 



H o<u 


! i.T 


l a 2Q 9 k 




h i<U 


i h iu 


! h 21,k 


! h 3U 


H 20,* 




r 

j_ H 22,* 


j_ n 32,A; 


H 30,A: 


! h 3i,* 


! H 32,* 


! H 33,* 



A l,0 T k 

8 i,o£* 



B 
B 
B 



B 



9A 

3,*_ 



[32] 



where H 0 o is a 3-by-3 matrix with the following form: 
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H 



00,/fc ~ Z 
9=1 



v +A * A Ll 



A L1 A * +A M 2 [C" 



'2 J -1 A 
,9 J A 



L2 A *y 



where Hi 0 is a l-by-3 row vector with the following form: 



h l0,* = 2 

9=1 



A, -ci T 



C f Ll -kiT A L1 A * ~ f L2 ' [ C v,? ]"' A L 2 A A: J 



where H 2 o is an N-by-3 matrix with the following form: 



H 20,* = 2 



9=1 



teN* +A L2 A L1 [ C nj]~ lA jt 
, " A L1 [ C " " A L2 A Il k^ 1 ^ 



kj 



where H 30 is a 2N-by-3 matrix with the following form: 



H 



3 

z 

9=1 



where hn is a 1-by-l matrix with the following form: 



41,* = £ 
9=1 



1 + c 2 -! 7 " 



L2y 



where h 2 i is an N-by-3 matrix with the following form: 



h 21,* = I 

9=1 



, + A £i J" 1 • f L1 + A L2 A£ 2 [cja J" 1 . f] 
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where h 3 i is a 2Af-by-l column vector with the following form: 



3 

I 
9=1 



LI 



9=1 



L2 



where H 2 2 is an N-by-N matrix with the following form: 



H 



22,k ~ ^ 
9=1 



[ C n}qf +A L2 A L1 [ C njT A L 2 A L1 
+ A L1 [c^V^ll + A L2 A L1 kir A L2 A L1; 



+ [»i,o]" 



where H 32 is a 2N-by-N matrix with the following form: 

f \ 



H 



32,* 



I -[ c "^' A i 

9 =A 



(-[ C v^ 1a L2 A Ll] 



3 

s 

9=1 



10 



where H 33 is a 2N-by-2N matrix with the following form: 



H 



33,/t 



i tear 


\ °NxN 






°NxN 


; 9=1 



+ D 



k-\ 



15 



where B 0 ,k is a 3-by-l column vector with the forms: 
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B 



0,k ~ Z 

H 



where Bi )k is a single value with the forms: 



10 



where B 2 ,i< is an N-by-1 column vector with the forms: 



*2,k = Z 
9 =1 



-A£l [cJi^A^p;" -A L2 A£? [c^^A^p^^ 



and where B 3k is a 2iV-by-l column vector with the form: 



B 3,* = 



7=1 \ J 



3 

z 

9=1 



(kiN* 



OP* 



+ D Jt-iAi, 0 N^r 



The matrix H A of form [32] is symmetric, and can be decomposed by a Cholesky 
1 5 factorization process into the form: 



H* = L* L* , where L* = 



L 00,A 


i 0 


i o 


! o 


'lO,* 




! o 


0 


L 20,* 


r 


| L 22,A: 


: ° 


L 30,A 


! *31,A: 


! L 32,A: 


! L 33,A: 



[33] 



-53 - 



10 



15 



Then, the unknown vector Y = 



8 1,0 1 ^ 

L A i,oN*j 



can be generated by a conventional forward 



and backward substitution processes by generating an intermediate matrix Y from the 
form L* Y = B k , and then generating matrix Y* from L* r Y* = Y , where 
B, 



B* = 



yA 
?A 

3 f k 



This generates an updated set of floating ambiguities and updated rover position. Next, 
matrix D* for the next iteration is generated as follows: D* = L33,* L33 3 * . The next 
iteration is then started by generating a new matrix H based on another epoch of data, 
and thereafter reiterating the above steps. While the epochs of data are generally 
processed in sequential time order, that is not a requirement of the present invention. In 
a post-processing situation, the epochs may be processed in any order. It can be shown 
that the matrix D* = L 33 ^ L33 7 * is mathematically equivalent to the following form: 



% = H 33,* -( H 30,£ ! h 31,* \ H 32,k) 



H 00,/fc 


! I.T 

| h io,* 


! H T 1 


-1 


a 30,k 


h l<U 


i V* 


1 1.T 

! h 2U 




"31,* 


H 2<U 


| h 2U 


! H 22,A j 




^ H 32,* 



[34] 



Matrix D* is the matrix of the second partial derivatives of the cost function 

F I, ($ x * > A 1,0 T * , 8 1>0 Ijt » A l,o N * ) partially minimized with respect to the variables 

8X *» A 1,0 7 *» 8 1,0 I * : 



D t =■ 



9[A 1<0 N^] 

where F 2 (A 1)0 N* ) = min F z (6X k , A lfi T k , 8 1)0 I* , A^qN* ) . 

8X *> Al.O 7 *. 6 !,!) 1 * 



[35] 
[36] 



The above process of generating the values of 8X k , Aj >0 r k , o 10 I^., and 
A 1,0 N * which minimizes the cost function F 2 (8X k , A 10 r A , 8 1)0 I ^ , A 1>0 N£ ) 
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generates an estimate for the position of Rover as Xq = (X 0 + 8X^ ), which has an 
accuracy which is acceptable for many applications. Higher accuracy generally results if 
we generate a set of fixed-integer ambiguities A 10 N^ from the floating ambiguities 
Al,oN£ , and generate an estimate of the rover position from the fixed-integer 
5 ambiguities. Examples of this are described in the next section. 

The means for performing then above steps in rover 100 are provided by main 
processor 110 under the direction of instruction sets stored in memory 112, with the 
various computed data being stored in data memory 114. 

10 Fixed-Integer Ambiguity Resolution 

Having at hand the estimation of the floating ambiguities A^qN^ , we describe 
an inventive method of generating the fixed-integer ambiguities 



' A l,oN L1 ^ 
l*l,oN L2j * J 



, which may be used alone or in combination with the processes 



described above. The components of the vector A 10 N L1 k have the same fractional part 
15 Ay/ Ll , and the components of vector A 10 N L2 k have the same fractional part At// 12 . 

Having at hand the matrix and the vector A 1?0 N^ , the embodiment forms the 
following cost function form: 

F(A W N*) = ±(A M N A -A l90 N^) r D^(A l50 N^ -A 1>0 N*) [37] 
The embodiment generates the values of A^qN^ which minimize the value of 
20 F(Ai^N k ) subject to the conditions: 

Ai,oN L1Jt eZ w * + Ai// Ll -I for some A^ L1 ,and 

A 1>0 N L2Jt g Z" k +Ay/ L1 • 1 for some Ay/ 12 , 

where Z nk denotes the domain of integer valued vectors (vectors having only integer 
components). To generate the values of A^qN^ , two reference satellites s Ll and s 12 
25 are selected for the first and the second frequency bands. A between-satellite double 
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difference matrix S is then generated, and will be applied to matrix A 1?0 N^ . Matrix E 
comprises a 2N-by-2N identity matrix, but with the columns associated with satellites 

^ j 1 r <-y 

s and s being modified, as shown below: 



2 = 















i 




1 


0 


-1 ... 


0 


0 


1 


-1 ... 


0 


0 


0 


1 ••• 


0 


0 


0 


-1 ... 


1 


0 



0 




V ! J 

The columns are modified by substituting a value of -1 for each column element that 
would normally be zero in the identity matrix. The positions of columns of matrix £ 
which are associated with satellites S 1 and £ 2 correspond to the row positions of these 
satellites in the vector A 1>0 N^ . When matrix S is multiplied onto vector A 1>0 N^ , the 

floating ambiguities corresponding to satellites S Ll and s L1 remain unchanged, but the 

floating ambiguity associated with s Ll is subtracted from the other floating ambiguities 

in the LI -band, and the floating ambiguity associated with s L2 is subtracted from the 
other floating ambiguities in the L2-band. Next, a permutation matrix II is generated, 
and is applied to matrix E to generated a matrix product II Permutation matrix is 
constructed to move the S x and £ 2 columns of matrix E to the first and second column 
positions of the matrix product II E. The construction of permutation matrices is well 
known to the field of mathematics. A permutation matrix satisfies the following 
relationships: II T -II = ITII T = I. Next, a change of variables for the cost function F is 



-56- 



undertaken as follows: 

Mft-n-£.A l9D Nt, 
M k =n.E.A 1?0 N*, 



[38] 



and the cost function can be formed as: 

F(M k ) = I (M 4 - M k f G k (Mjfc - M* ) . [39] 

The first two components of are real valued components (floating point numbers) 
and the remaining 2N-2 components are integers. Let us divide the matrix G k into 



appropriate blocks in accordance with division of the vectors M k = 



l,k 



and 



M t = 



M 

Mi 



l,k 



2,kJ 



G ll,£ ! G 21,A: 
1- 



G 21,£ ! G 



22,*; 



where G n k is a 2-by-2 matrix, is a 2-by-(2AT-2) matrix, and where G 2 2,A: is a 

(2N-2)-by-(2N-2) matrix. This partitioning enables us to use a two step process of 
minimizing F(M k ) = ± (M k -M k ) T G k (M k - ) by first generating values of 

M l k that reduce the value of F(M k ) , and thereafter generating values of M 2fk that 
further reduce the value of F(M k ) . We notate this two step process as: 

}- 



mm 



{i(M^M^) r G^(M^-M^) 



mm min<! 

M 2( a Mi,/t 



±(M hk -M M ) T G 1U (M U -M 1>t ) 
+ (M 2>/t -M 2)/t ) T G 2U (M u -M U ) 
+ l(M 2>/t -M 2>it ) T G 22yk (M 2>k -M 2tk ) 



[40] 



The first step (inner minimization) of generating Mj ^ to reduce F(M^.) is performed 
in the space of floating point numbers (real valued space), and generates M 1)t in a form 
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equivalent to: 



M U =M U +[G lu ]" 1 [G 2U f(M u -M 2fk ) 



[41] 



the second step (outer minimization) of generating M 2> £ is performed by substituting 

the form [41] of Mj k into form [40] to generate the following modified version thereof: 

5 {(M 2) * -M 2 ,*) r [G 22) * -G 2l9k [G n9k ]- l [G 21ik ] T ](M 2ik -M 2 ^ k ) [42] 

Then, a search of an integer subspace for values of M 2 ^ is performed to find a set of 
integers which provides the smallest value of the modified form [42], The integer 
subspace can be relatively small, and centered about the floating point values of M 2 k . 

After the integer search of form [42] is performed, the result M 2> £ is substituted into 

10 form [41] to generate a revised vector M l k . Then, the change of variables operation 

(form [38]) is reversed as follows to generate the fixed-integer ambiguities A lj0 N^ as: 

A h0 N k =H l n T M k [43] 

A 1>0 N£ may then be separated into its components Aj^N^ 1 and A^qN^ 2 . 

As an option, an ambiguity resolution validation procedure may be performed to 
15 check for consistency in the ambiguity resolution. This procedure is conventional, and 
the reader is referred to prior art literature for a description thereof. 

Having resolved the fixed-integer ambiguities, we can new generate a further 
refined estimate of the other solvable unknowns: 5X^, A 10 r^, d 10 l k . A third cost 
function is formed as follows: 

20 F 3 (dX k , A li0 r k , & h0 I k ) = 



(a^X* -cA 1)0 r* 1 +5 lj0 I* -A q , 0 l>; L1 f [c^]" 1 > ^ 
, (a*SX*. -cA 1j0 t* 1 +8 lj0 Iyt -A qj0 b* L1 ) 



>{ (a*SX* -cA lfi T k -1 + Al 2 A-2 8 1)0 I* -A q , 0 b; L2 ) 



+ ^*«A -cA h0 r k -1 + A*2 Al] 5 h0 I k -A^b^f^SjJ" 1 x ^ 
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where the five terms of the form are the same as the second through sixth terms of form 
5 [31], except that floating ambiguities have been replaced by the fixed-integer 

ambiguities. The same estimation procedures used on form [31] may be applied above to 
form [44], except that only one iteration is needed, and matrix D is not generated. 

The means for performing then above steps in rover 1 10 are provided by main 
processor 110 under the direction of instruction sets stored in memory 112, with the 
10 various computed data being stored in data memory 114. 

If lower accuracy can be tolerated, we may apply the process on form [44] using 
the floating point ambiguities rather than the fixed-integer ambiguities. In this case, the 
above fixed-integer ambiguity resolution process may be omitted. 

15 Between Base Station Processing, Part II 

In the above embodiments, the data of forms [6] and/or [7] were provided to the 
rover station. However, because several of the rover-to-base station processes described 
above are more efficient than those in the prior art, it is possible for the rover itself to 

20 undertake the task of generating some or all of forms [6] and [7] at the rover location, in 
real-time, from phase and pseudo-range measurements conveyed from the base stations 
to the rover. This information may be conveyed by radio-signals from the base stations 
to the rover, as described above. One may also implement a system whereby the base 
stations convey their information to a relay station by cable (such as the internet), with 

25 the relay station positioned within a few kilometers of the rover station. The relay 
station then relays the base station data to the rover by radio-signals. 

In these embodiments, the rover receives the following data from each base 

station "r" (r = 1, 2, 3) for several epochs "k": pQ 5 , pffi , fQ s , and </>^ s . 
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The locations of the base stations are also obtained: Xi, X 2 , and X3. As indicated 
previously, the locations of the satellites are highly predictable and can readily be 
determined by the rover with its clock and correction data from the almanac data 

transmitted by the satellites. From this, the rover generates the computed ranges R s . 

~ *-' r,K 

from each rover "r" to each satellite "s" of a group of satellites being tracked. As before, 
the troposphere delay terms are estimated from the Goad-Goodman model. 

With this information, the residuals (difference quantities) A qr b^ ,s , 

A g,r b k 2 ' S ' A q,rPk X ' S > 311(1 A^,./?""* of forms [4A-4D] are generated for each 
baseline (q, r), where (q, r) has the following pairings (B2, Bl), (B3, Bl), and (B2, B3). 
For each baseline, the solvable unknowns in forms [5A-5D] are A q r r k , A qr I k , and 

Aq r N^ . Values are estimated from the residuals in a manner similar to that described 

about for the primary baseline between the Rover and first base station. As an example, 
a cost function F(*) may be formed as follows: 

F ( A q,r r k> &q,r l k> A q,r N A:) 

= |(A q , r N >fc _ 1 -A q>r N Jfc ) r D ft _ 1 (A q , r N it _ 1 -A q>r N*) 



+ i 



( 



{-cA q ,r k l+A^I* -A q , r b^f lp™ >r Y x 
< (-cA^-l+Aq^-AqX' 1 ) 

' (-cA q>r r k i + Al 2 A" 2 ■ A w I t -A q4 .bJ 2 f [c^]" 1 x ' 
, (- cA 1,rn ■ I + A L A-2 • A q ,I A - A q;r b^ 2 ) 

(-f Ll A q>r r k -1 - A£\ .A q , r I, +A q>r N L1 -A q , r p^ J[c^ q>r \ l x ' 
, (-f Ll * g ,rn 'I -A£\ A q>r I A + A q , r N L1 -A q , r p^ ) j 

f L2 A q ,r k 1 -A L2 A" 2 • A q , r I* +A q;r N L2 -A q?rP £ 2 f [c^]" 1 x > 
L2a _ 7 * * -2 A T . a T?L2 a _L2 j 

[45] 



, (-f L2 A g , r r k 1 -A L2 A" 2 A q>r I A + A q , r N L2 -A q?r p£ 
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wherein the covariance matrices C^J r , C^ q r , c\J q r , C , a q r are generated in a 

similar manner as the covariance matrices for the primary baseline. 

The floating ambiguities may be estimated in an iterative manner as described 
above. We start at the initial epoch £=0 with the weighting matrix D 0 set to the zero 
matrix, and an initial guess of floating ambiguities A q9r N 0 equal to zero. The first term 
of form [31] evaluates to zero for this initial epoch. We then generate a set of values for 
Aq 9r Tk, A q>r l£, A q5r N^ at a first epoch k=l which moves the value of F towards 

zero. This generates an initial estimate Aq^Nj for the floating ambiguities. A 
weighting matrix Di for the ambiguities is then generated, and a new set of values is 
generated for A^ r r^, A q?r I^, A qjr N^ at the next epoch k~2 which moves value of F 

(k=2) towards zero. As a result, subsequent estimates A q?r N 2 is generated. This 
iteration process continues, with Aq^N^ generally improving in accuracy as the 

iterations progress. The corresponding H-Y = B form for this process (similar to form 
[32]) is as follows: 





| h 21,k 


l h 31,k 


h 21,k 


! H 22,k 


i¥¥ T 


h 31,k 


! H 32,/t 


i 

! H 33,* 





A q,r T k 












?A*_ 




_Aq,rN*_ 




B 3,k 



[46] 



where hn is a 1-by-l matrix with the following form: 

K + tll-fcwrY f Ll + f L2 •[c",r] _1 -f L2> 



41,* 



where h 2 i is an N-by-3 matrix with the following form: 



l 2l,k 



-c ^fwY 1 " e Al 2 A£? [c\%, r Y 1 
+ A£i [c^J 1 f L1 + A L2 All [c^,rV ' L 2> 
where I131 is a 2jV-by-l column vector with the following form: 
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10 



l 3l,k 



L2 



where H22 is an N-by-N matrix with the following form: 

H 22,* = knlrV + A L2 A M fe^L A L1 + A L1 fcSi,]"^ 



+ A ¥ • A 



-2 UL2 



2 t -1 



—2 — - 
L2 ^Ll L^v^r J A L2 A L1 + t a q,r] 



r - n-2 



where H 32 is a 2N-by-N matrix with the following form: 



H 



32,* 



-[ C vi,rf l A L 



{-fe^U A L 2 l) 



where H 33 is a 2N-by-2N matrix with the following form: 



H 



33,* 



fc L1 V 




°NxN 


®NxN 




[c L2 V 



D 



k~\ 



where Bi, k is a single value with the forms: 



1 5 where B 2 ,k is an JV-by-1 column vector with the forms: 



B 2,k ~ 



[C^r]" 1 A q , r b™ + Ah All HXrV^Y 
-A£j kV.^A^p^ -A L2 All teXrV KrPk 2 , 
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and where B 3)k is a 27V-by-l column vector with the form: 

f \ 
-l 



,L1 H 



fc L2 Y 



A q,rP^ 



+ D *-l A q,i-N*_r 



The matrix H* of form [46] is symmetric, and can be decomposed by a Cholesky 
factorization process into the form: 







0 0 




— L* L* r , where = 




| L 22,/fc 


0 


[47] 




l 31,k 


! L 32,/fc 


L 33,/t 





Then, the unknown vector Y = 



A q,r T k 

A q,rl* 
Aq,rN* 



can be generated by a conventional forward 



and backward substitution processes by generating an intermediate matrix Y from the 
form L* Y = B k , and then generating matrix Y* from L, k T Y* = Y , where 



8* = 



2,* 



B 



3, A: 



This generates an updated set of floating ambiguities and updated rover position. Next, 
matrix D k for the next iteration is generated as follows: D* = L 33 ,* L33 7 * . The next 
iteration is then started by generating a new matrix H based on another epoch of data, 
and thereafter reiterating the above steps. While the epochs of data are generally 
processed in sequential time order, that is not a requirement of the present invention. In 
a post-processing situation, the epochs may be processed in any order. It can be shown 
that the matrix D* = L 33 ,* L33 7 * is mathematically equivalent to the following form: 

r 1 V 1 



D * = H 33,* - ( h 31,* ! H 32,* ) 



h 1 h A 

^ h 21,* ! H 22,A J 



H 32,* J 



[48] 
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Matrix Dfc is the matrix of the second partial derivatives of the cost function 
F (&q,r T k > Aq,!- 1 * > A q,r N & ) partially minimized with respect to the variables 

A q,r T k> &q,r l k 

*2 



[49] 
[50] 



where F(A q>r N A ) = minF (A 9>r r*, A^I*, A q>r N^) . 

The corresponding fixed-integer ambiguities A q>r N A may be generated from the 
floating ambiguities A q>r N£ by the same process described above with reference to 

forms [37] - [43] that is used to generate the fixed-integer ambiguities associated with 
the primary base line between the rover (R) and first base station (Bl). 

Having resolved the fixed-integer ambiguities, we can new generate a further 
refined estimate of the other solvable unknowns: A q r z k and A qr I^.. A third cost 
function is formed as follows: 
F 4 ( A 9) r^/t . Aq,!- 1 * » A q,r N A ) 

J (-cA qjr T k -i + A qtr I k -A q>r b™Y [ C nt,,rY X ' 
{ (-^ ;r r A -l + A q>r I jt -A q , r bJ' 1 ) j 

A (-cA ff ,r* 1 + Al 2 AJ" 2 A q ,I k -A q>r bJ 2 f [ctj, r ]" ! * 



1 + A 



-2 



L2 A L1 A q,r l k 



L2 



) 



-1 - A£j .A V I 4 +A q , r N L1 -A q>r p, L1 f [c^]" 
, (-/"a,,t 4 -1 - A"\ .A V I 4 + A q , r N L1 -A^pJ* ) 



•1 ^ 

x 



(-/"A,,r* 1 - A L2 A" 2 A^I* + A q , r N L2 -A q , r p^ f [ C ^ r \ X x ^ 
, (-/"A,, r r, 1 - A L2 A£ 2 • A q , r I A + A q , r N L2 - A q , rP ^ 2 ) 

+ ^ q>r I 4 r [o q>r ]- 2 A q>r I A [51] 

where the five terms of the form are the same as the second through sixth terms of form 
[45], except that floating ambiguities have been replaced by the fixed-integer 
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ambiguities. The same estimation procedures used on form [45] may be applied above to 
form [51], except that only one iteration is needed, and matrix D is not generated or used. 
As a result, estimates for A^ r r^ and A q r l£ are generated. However, we prefer to 

perform some consistency checks on these estimates before providing to the process that 
5 operates on the primary baseline between the rover and the first base station. Thus, we 
will denote these estimates with hat symbols as follows: A g r f k and A qjI .i^ . After the 

estimates associated with the three base line (q,r) = (B2, Bl), (B3, Bl), (B3, Bl) have 
been generated for the k-th epoch, we have the following data: 

A 2,l^t> A 2,li*> A 2,lN*> 

10 A 3>1 f*, A^I*, A 34 N^,and 

A 3,2^, A 3,2**> A 3,2Na;. 
We first perform an ambiguity resolution closure check. In this check, the following 
relationships should hold: 

A 2,lN* + A 3 , 2 N* + A lj3 N* =0, [52A] 
1 5 which is equivalent to : 

A 24 N* + A 3 , 2 N* - A 34 N* = 0 . [52B] 

If the above relationships are not satisfied, the ambiguities have been resolved 
incorrectly for at least one base line, and the estimations of fixed ambiguity should be 
neglected. To address an incorrect resolution, new data may be taken, various subsets of 
20 data may be processes to generate sets of ambiguities which satisfy the above 

relationships. After all three between base receivers ambiguity have been resolved and 
the above relationships are satisfied, the ambiguity vectors A 2jl N£ , A 3?1 N^ , and 

A 3,2^£ are considered to be correctly fixed. 

The time delays are similarly checked and corrected by taking new data or 
25 searching existing data sets. In this check, the following relationships should hold: 
A 2,l^£ + A 3,2^£ + A l,3^A: = 0 , to within a tolerance value of ±G\ , [53A] 
which is equivalent to: 

A 2,l*> + A 3,2*> *~ A 3,l^Jfc = 0 , to within a tolerance value of ±€\, [53B] 

where 6 1 is a tolerance level which is close to zero. In general, 6 \ depends upon the 

-65- 



distances between baselines and the desired degree of accuracy for the system. If the 
above relationships hold, then we validate these delays (i.e., A 21 T£ = A 21 f k , 

A 3,l^ = A 3,l f k > and A 3 2 t^ = A 3>2 f k ). 

Next, we perform a consistency check on the ionosphere delays. In the following 
discussion, we the subscripts on the ionosphere delay differentials have sometimes been 
exchanged, but this is of no substantive consequence since Aj^I^ = -A 2)1 l£ , 

ft * A * A 

A 2,3^A: = ~A 32 l£, Aj^I^ = -A 31 l£. The between base receivers ionosphere 
estimations should satisfy the relationship 

* /\ /\ 

^l,! 1 * + ^2,3!* + ^3^1^ = 0 , to within a tolerance value of ±€ 2 . [54] 
However, measurements noise typically prevents these relationships from being satisfied 
to an acceptable tolerance level ±6 2 . To better satisfy [54], the following quadratic 
function is minimized to obtain new estimations A lj2 ^k'^2^k>^3,i^k '• 
F I ( A l,2 1* . A 2,3 h . *3,lh ) 

= \^i,ih -A 1>2 4) r [c lj2 r 1 (A u T t -A 1>2 i*) 

+ \(*2jh -*2jh) T [c 2 ,3r\±2jh -A 2 , 3 i/t) [55] 

— > _ min _ , 
provided the condition of: 

*iah + ^2,3^ + ^3,1^ = 0 [56] 

is satisfied. The ionosphere estimations covariance matrices C 12 »C 23 ,C 35 j are 

estimated in the previous process by conventional methods. The inventors have found 
that the following form minimizes form [55] subject to condition [56]: 

±i 9 ih =*i 9 iik -Ci, 2 [c u +c 2>3 +c 34 r 1 (A 1?2 i* +A 2j3 i* + a 34 i*) 

*2jh =*2jh -C 23 [C lj2 + C 2>3 +C 34 r 1 (A l92 I^ +A 23 I^ +A 34 I*) [57] 

= A 3a i t -c 34 [c w +c 2j3 +c 34 r 1 (Ai, 2 i^ +a 23 i^ + A 3>1 i*) 

Before applying form [57], the vectors &\^k>&2^k^3j\k are preferably smoothed 

using Kalman filtering scheme with a dynamic model based the above described the 
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Gauss-Markov time model (forms [GM1] and [GM2] above). 

The results of these process may be provided to the process which generates the 
ionosphere delay differentials for the baselines associated with the rover stations, 
specifically the differentials generated according to forms [26], [29], and [30]. For form 

[26], which generates A^qI^ , the results are provided as follows: A 1?2 l£ = A 12 l£ 

Ai^I* = ^l^k • For generating the ionosphere delay differentials according to form 

[29] we set A 21 l£ = ~Ai,2 l k • For generating the ionosphere delay differentials 

according to form [30] we set A 3)1 I^ = A 3)1 T^ . 

Thus, to summarize a preferred embodiment of the present invention, at each 
epoch the following occurs for the three baselines between the base stations: 

• Update of the floating ambiguity estimates by way of forms [45]-[50] , 

• Generate fixed-integer ambiguities (with a process similar to that illustrated by 
forms [37] - [43]); 

• Generate estimates of A qr r k and A qr I k by way of form [51] with a process 
similar to that illustrated by forms [45]-[50] 

• Check consistency of the ambiguities (\ 2 ,l^k + A 3>2 N yt + Aj^N^ = 0) and 
time offsets (A 2 ,l*k + A 3j2 t* + A^f* = 0); 

• Generate adjusted estimates A 1>2 1* , A 2 3 T^ , A 31 of the ionosphere delay 
differentials; 

• Provide these results to the process of estimating the rover's position using the 
primary and secondary baselines. 

The above base-station to base-station data can be generated by the rover, or by 
an external source, such as relay station. In the case of rover 100, the means for 
performing all of the above steps are provided by Between-Base Station processor 140 
under the direction of instruction sets stored in an instruction memory, with the various 
computed data being stored in a data memory. Moreover, while higher accuracy is 
obtained by generating the fixed-integer ambiguities, lower accuracy embodiments may 
just generate the floating ambiguities. 

While it is preferable to interpolate the ionosphere delays using two baselines 
between three base stations, it may be appreciated that some application may achieve 
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acceptable accuracy by only interpolating the ionosphere delays using one baseline 
between two base stations. Such an example may be a road project where the road is 
relatively straight, as shown in FIG. 4. 

Each of the above methods of generating the base station data and estimating the 
5 coordinates of the rover is preferably implemented by a data processing system, such as a 
microcomputer, operating under the direction of a set of instructions stored in computer- 
readable medium, such as ROM, RAM, magnetic tape, magnetic disk, etc. All the 
methods may be implements on one data processor, or they may be divided among two 
or data processors. 

10 

Computer Program Products 

It may be appreciated that each of the above methods may comprise the form of a 
computer program to be installed in a computer for controlling the computer to perform 
the process for estimating the location of a rover station (R) with the use of at least a first 

15 base station (Bl) and a second base station (B2), with the process comprising the various 
steps of the method. 

In addition, each of the above methods of generating the base station data and 
estimating the coordinates of the rover may be implemented by a respective computer 
program product which directs a data processing system, such as a microcomputer, to 

20 perform the steps of the methods. Each computer program product comprises a 

computer-readable memory, such for example as ROM, RAM, magnetic tape, magnetic 
disk, etc., and a plurality of sets of instructions embodied on the computer-readable 
medium, each set directing the data processing system to execute a respective step of the 
method being implemented. FIG. 7 shows an exemplary comprehensive listing of 

25 instructions sets for implementing the above method. Each of the above methods is 

achieved by selecting the corresponding groups of instruction sets (as apparent from the 
above discussions). Instruction set #18 is common to all of the methods, and is modified 
to omit the data which is not used by the particular method. 

While the present invention has been particularly described with respect to the 

30 illustrated embodiments, it will be appreciated that various alterations, modifications and 

adaptations may be made based on the present disclosure, and are intended to be within 

the scope of the present invention. While the invention has been described in connection 

with what is presently considered to be the most practical and preferred embodiments, it 
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is to be understood that the present invention is not limited to the disclosed embodiments 
but, on the contrary, is intended to cover various modifications and equivalent 
arrangements included within the scope of the appended claims. 
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